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INTRODUCTION 


This  Annual  Scientific  Report  covers  work  performed  under 
Contract  No.  F19628-75-C-0050;  previous  related  contracts  were 
F19628-72-C-0121,  AF19(628)-5187,  and  F19628-69-C-0119. 

The  work  comprises  experimental  and  theoretical  determinations 
of  cross  sections  and  reaction  rates  for  the  mutual  neutralization  of 
positive  and  negative  ions  of  the  type 

a'*’  +  b"  -  A  +  B  . 

These  reactions  play  an  important  role  in  determining  the  ion  density 
and  composition  in  the  D-region  of  the  earth's  ionosphere.  This 
reaction,  along  with  photodetachment,  photodissociation,  and  positive 
ion-electron  recombination,  terminates  the  chain  of  events  following 
the  initial  photoionization  of  the  earth's  atmosphere  between  about 
60  and  95  km.  The  reaction  is  also  important  at  lower  altitudes  when 
the  atmosphere  is  disturbed,  such  as  in  the  ionization  surrounding  a 
nuclear  explosion. 

Most  of  the  neutralization  reactions  between  simple  atmospheric 

ions  were  studied  quite  successfully  under  previous  contracts. 

Experimental  measurements  were  supplemented  by  theoretical  calculations 

(without  contract  support)  which  led  to  a  fairly  good  understanding 

of  the  reaction  mechanisms.  More  recently,  however,  the  ionospheric 

problems  hi've  centered  on  reactions  between  water-clustered  ions  of  the 

type  NO^.  (H-0)  ,  H^.  (H„0)  ,  C0_  .  (H.O)  ,  and  NO,  .  (H,0)  ,  where  n  is 
2n  2nj  2n  IZn 

usually  between  1  and  about  4.  These  ions  are  formed  slowly  by 
three-body  association  reactions  and  dominate  the  quiescent  night-time 
D  region,  and  probably  prevail  both  day  and  night  at  altitudes  below 
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50  km.  Because  of  the  slow  rates  of  their  formation,  and  the  weakness 
of  the  cluster  bond  (about  1/2  eV),  these  ions  are  very  difficult  to 
form  in  beams  such  as  are  required  by  our  merging-beam  apparatus. 

Because  we  anticipated  experimental  difficulties  and  delays,  and  because 
there  were  some  important  physical  differences  between  the  neutralization 
reactions  with  clustered  ions  and  with  the  more  simple  ions,  we 
allocated  a  significant  ^unount  of  our  contractual  effort  toward  theory. 
This  effort  has  yielded  important  insights  into  the  reaction  mechanisms, 
and  is  approaching  the  capability  of  actual  cross  section  and  rate 
calculations. 

The  current  state  of  the  work  is  described  below.  In  addition, 

we  are  submitting  three  journal  reprints  as  part  of  this  report.  One^ 

2-3 

relates  directly  to  the  current  work,  and  the  others  concern  work 
performed  under  previous  contracts,  but  were  published  during  this 
contract  period. 
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EXPERIMENTAL  RESEARCH 


The  experimental  effort  has  been  directed  primarily  toward 
improvements  to  and  testing  of  the  apparatus.  Several  ion  source 
configurations  that  showed  promise  for  cluster-ion  production  have  been 
tried.  Toward  the  end  of  the  last  contract,  a  low  temperature  hollow 
cathode  g\ow  discharge  was  studied.  It  produced  hydrated  ion  currents 
of  small  but  promising  magnitudes.  Under  this  contract  it  was  found 
that  the  energy  spreads  in  the  negative  ion  beams  from  this  source 
were  prohibitively  large  and  prevented  the  achievement  of  sufficient 
quality  in  the  beams.  Analysis  of  the  ion  optics  in  the  apparatus  also 
showed  several  weaknesses  that  decreased  the  signal/noise  characteristics. 

Improvements  were  undertaken  under  this  contract  along  three  basic 
lines.  Firstly,  to  produce  hydrated  ion  beams  of  reasonably  high 
intensity  and  narrow  energy  spread,  we  built  and  tested  hollow  cathode 
duoplasmatron  sources.  The  low  currents  of  hydrated  ion  beams  produced 
even  by  these  highly  efficient  sources  required  an  improvement  in  the 
sensitivity  of  the  apparatus.  For  this  second  goal,  we  made  changes 
in  the  ior.  optical  system  to  improve  the  ion  betun  quality,  the 
exactness  with  which  the  ion  beams  are  merged  and  demerged,  and  the 
sensitivity  with  which  the  neutral  beam  can  be  detected.  Thirdly, 
the  vacuum  system  was  reorganized  and  poor  components  were  repaired  or 
replaced  in  order  to  improve  the  vacuum  attained  as  well  as  to  increase 
ease  of  operation  and  to  decrease  the  chance  of  a  catastrophic  vacuum 
failure  through  the  use  of  various  automatic  vacuum  safety  devices. 
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THEORETICAL  RESEARCH 


lon-Dtpole  Ion  Scattering  Theory 

During  Che  first  year  of  this  contract  we  have  continued  the 

theoretical  study  of  Che  mechanism  of  recombination  of  hydrated  or 

clustered  ions  in  Che  ionosphere  chat  was  begun  under  the  preceding 
4 

contract  period.  To  date  we  have  focused  our  attention  on  the  inter¬ 
action  of  an  atomic  ion  with  a  rigid  dipolar  ion  as  a  model  problem 
for  developing  and  testing  approximation  techniques  Co  be  applied 
to  Che  more  complex  situation  involving  clustered  ions. 

In  the  initial  phase  of  this  effort  we  studied  Che  ion-dipolar 
ion  problem  in  two  dimensions.  Using  Che  efficient  and  rapid  comput¬ 
ing  scheme  that  we  developed  earlier  for  integrating  Che  fixed  para¬ 
bolic  orbit  classical  equations  of  motion,  we  have  carried  out  a 
systematic  investigation  of  the  trajectories  to  determine  the 
preferential  conditions  leading  to  energy  transfer  and  orbital  Crapping. 
We  have  found  that  the  trajectories  can  be  characterized  by  the  number 
of  times  Che  rotor  reverses  direction  as  the  ion  goes  by,  this  number 
depending  on  such  orbit  parameters  as  Che  disCtince  of  closest  approach, 
the  relative  phase  of  the  rotor's  motion,  and  the  initial  rotor  energy. 
As  the  conditions  for  rotor  reversal  are  approached,  large  energy 
transfer  becomes  possible.  The  capture  cros*^  section  depends  on  the 
relative  amount  of  phase  space  associated  with  the  transitions  giving 
large  energy  transfer.  Figure  1  shows  Che  typical  results  under 
conditions  of  strong  interaction  between  the  ion  and  the  charged  dipole. 
Here,  the  rotor's  reduced  angular  velocity  (omega)  has  been  plotted 
versus  the  angular  position  of  the  ion  4)hi  ■  -  n  initially  and  ■  0 
at  Che  turning  point  of  the  trajectory)  for  a  wide  range  of  initial 
rotor  phases,  all  other  initial  conditions  remaining  constant. 
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We  have  carried  out  a  systematic  comparison  of  these  fixed  orbit 
numerical  calculation  results  with  those  obtained  fr^^m  a  first  order 
perturbation  solution  of  the  equations  of  motion  in  which  the  rotational 
energy  transfer  is  expressible  in  closed  form  in  terms  of  Airy  functions 
and  their  derivatives.  These  comparisons  suggest  that  the  simple 
perturbation  approach  is  a  reasonably  good  first  approximation  to  the 
maximum  energy  transfer  available  under  given  conditions.  These  results 
x;ere  reported  by  Dr.  Felix  T.  Smith  in  a  paper  contributed  to  the 
meeting  of  the  Division  of  Electron  and  Atomic  Physics  of  the  American 
Physical  Society  in  Chicago,  December  2-4,  1974. 

Following  the  development  of  satisfactory  techniques  for  treating 
the  two-dimensional  ion-dipolar  ion  problem,  we  turned  our  attention 
to  extending  the  number  of  degrees  of  freedom  that  we  could  handle. 

To  begin  with,  we  concentrated  on  the  extension  from  motion  in  a  plane 
to  motion  in  three-dimensional  space,  which,  for  our  problem,  brings 
in  two  additional  degrees  of  freedom.  We  initially  based  this 
extension  on  the  analogous  work  of  R.  J.  Cross,  and  in  this  context  we 
reviewed  and  extended  his  work  in  connection  with  the  interaction  of 
an  ion  with  a  neutral  dipole.^  His  perturbation  theory  approach  is 
closely  analogous  to  the  one  we  have  ^een  using  except  that  the  problem 
of  the  neutral  dipole  involves  trajectories  that  are  straight  lines 
instead  of  parabolic,  with  the  result  that  the  Bessel  functions  and 
appear  instead  of  the  Airy  functions  that  are  appropriate  to  the 
coulomb  problem.  To  determine  the  capture  probability,  which  is  a 
function  of  the  initial  relative  kinetic  energy,  initial  dipole  rotational 
energy,  and  the  impact  parameter,  it  is  necessary  to  carry  out  an 
averaging  procedure  over  the  three  angles  (specifying  the  orientation 
of  the  rotor's  position  and  angular  velocity  at  some  reference  point) 
appearing  in  the  expression  for  the  rotational  energy  transfer. 
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Unfortunately,  we  have  found  that  the  necessary  averaging  procedure, 
though  possible  in  principle,  would  require  an  excessive  computational 
effort. 

To  avoid  this  difficulty  we  reexamined  the  problem,  making 

maximum  use  of  a  representation  in  angular  momentum  rather  than  angular 

coordinates  so  as  to  take  advantage  of  angular  momentum  conservation 

laws.  In  this  context  we  have  been  able  to  make  use  of  the  work  of 
6,7 

Miller,  whose  approach  leads  to  a  considerable  simplification  of  the 
problem  and  allows  treatment  of  the  full  three-dimensional  motion  with 
an  apparatus  that  is  hardly  more  complicated  than  is  needed  for  the 
two-dimensional  problem.  We  made  use  of  the  first-order  classical 
perturbation  treatment  that  led  to  an  expression  in  terms  of  Airy 
functions  in  the  two-dimensional  case  and,  as  it  turned  out,  exactly 
the  same  integrals  suffice  for  the  new  three-dimensional  expressions. 
The  perturbation  theory,  when  applied  semiclasslcally,  gives  us  a 
convenient  closed-form  expression  for  the  full  scattering  matrix, 
and  Professor  Miller  assisted  us  in  extending  this,  using  a  uniform 
approximation  to  cover  quantum  effects  in  the  nonclassical  region  of 
the  motion.  The  only  limitations  of  this  procedure  apply  in  very 
close  encounters  when  the  interaction  forces  are  so  great  that  first 
order  perturbation  theory  ceases  to  be  appropriate.  We  feel  that  these 
will  be  for  very  small  classical  impact  parameters  and  ordinarily  will 
not  contribute  much  to  the  cross  sections  of  interest.  In  any  case, 
they  can  be  treated  through  the  use  of  classical  trajectory  calcula¬ 
tions  for  which  we  have  already  developed  good  computing  programs. 

A  paper  describing  our  semiclass ical  perturbation  technique  was 
published  in  the  October  20,  1975  issue  of  Physical  Review  Letters.^ 

A  copy  of  that  article  is  included  as  Appendix  A  of  this  report. 
Furthermore,  details  of  the  new  semiclassical  technique  and  some 
preliminary  results  on  the  cross  sections  for  capture  of  0  by  NO^  were 
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presented  by  Dr.  Felix  T.  Smith  in  a  contributed  paper  at  the 
IXth  International  Conference  on  the  Physics  of  Electronic  and  Atomic 
Collisions  held  at  the  University  of  Washington  in  Seattle,  Washington, 
July  24-30,  1975. 

Since  that  meeting,  considerable  effort  has  been  devoted  to 
improving  and  checking  the  computing  schemes  for  evaluating  the 
capture  cross  sections  starting  from  the  perturbation  semiclassical 
form  of  the  S-matrix.  Thus  we  are  now  in  a  position  to  generate  in  a 
systematic  fashion  further  capture  cross  sections  for  reactions  of  the 
type  A  +  BC^  where  the  molecular  ion  has  varying  dipole  moment  and 
varying  moment  of  inertia. 

Related  Work 

As  explained  in  the  reprint  of  our  Physical  Review  Letters  article 

(Appendix  A),  our  semiclassical  perturbation  approach  is  equally 

applicable  to  treating  the  interaction  of  an  electron  with  a  neutral 

dipole.  Under  the  sponsorship  of  other  agencies,  and  with  SRI  Internal 

research  funds,  we  have  devoted  considerable  effort  to  computing  the 

elastic  and  rotationally  Inelastic  cross  sections  for  electron  impact 

on  highly  polar  neutral  molecules  such  as  CsF  and  KI.  Our  interest  in 

these  two  systems  arose  because  these  cross  sections  have  been 

8  9 

measured  experimentally  by  Stern  and  his  students  ’  at  Columbia 

10 

University  and,  furthermore,  they  have  been  and  currently  are  the 
subject  of  study  at  other  laboratories  using  the  close-coupling  quantum 
mechanical  approach. 

We  have  now  used  the  semiclassical  S-matrix  to  compute  the  total 
scattering  cross  section  (elastic  and  rotationally  inelastic)  for  a 
test  case  in  the  scattering  of  an  electron  by  CsF  at  1  eV  with  the  polar 
molecule  Initially  in  the  rotational  state  with  j  =  41.  We  have  been 
able  to  compare  the  results  with  calculations  by  Norcross  and  Collins 
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at  JILA  (Boulder,  Colorado)  who  have  done  it  with  a  Born  approximation 
and,  for  smaller  quantum  numbers,  the  close-coupling  calculations. 

Our  results  deviate  significantly  from  the  close-coupling  results  only 
at  very  small  values  of  orbital  angular  momentum  quantum  numbers, 
namely  for  A  ^  2.  In  this  region  the  Bom  approximation  is  totally 
inadequate  and  it  does  not  converge  to  the  semiclass  leal  values  until 
about  X  =  7. 

The  ability  of  the  perturbation  semiclassical  approach  to  success¬ 
fully  treat  the  electron-polar  molecule  Interaction  problem  at  all  but 
the  smallest  quantum  numbers  demonstrates  the  great  flexibility  of 
our  method.  It  is  particularly  noteworthy  that  our  method,  in  contrast 
to  either  the  exact  quantum  mechanical  or  the  fully  classical  approach, 
enables  us  to  understand  both  qualitatively  and  quantitatively  the 
contributions  corresponding  to  essentially  classical  scattering,  to 
quantum  interference  effects,  and  to  tunneling  into  a  classically 
forbidden  region.  As  for  the  failure  of  the  semiclassical  method  at 
the  few  lowest  lying  orbital  angular  momentum  quantum  numbers,  we 
believe  that  there  is  considerable  scope  for  refining  our  approach  so 
as  to  reduce  the  discrepancy  between  the  close  coupling  and  semiclassical 
results  in  this  quantum  number  regime.  This  could  be  achieved  by 
(a)  removing  the  perturbation  approximation  and  instead  using  the 
numerical  solution  of  the  exact  classical  equations  of  motion  of 
the  system  for  computing  the  S-matrix  elements  and  (b)  using  a  more 
realistic  short  range  potential  by  supplementing  the  charge-dipole 
Interaction  term  with  higher  order  multipole  and  polarizability  terms. 

Returning  to  the  ion-ion  mutual  neutralization  problem,  we  feel 
that,  due  to  the  long  range  nature  of  the  coulombic  forces,  the  total 
orbital  capture  cross  sections  are  likely  to  be  dominated  by  contribu¬ 
tions  from  trajectories  characterized  by  relatively  large  distances  of 
closest  approach.  The  results  of  our  preliminary  calculations  for  the 
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capture  of  0  by  NO  lend  strong  support  to  this  view.  Hence  we 
expect  the  pert”rbatlon  semlclasslcal  approach  as  it  stands  to  be 
adequate  for  thi?  calculation  of  reliable  lon-lon  neutralization  cross 
sections.  The  considerations  relating  to  electron- polar  molecules  dis¬ 
cussed  in  the  previous  paragraph  are,  however,  relevant  to  the  question 
of  calculation  of  electron-ion  recombination  cross  sections  and  we 
anticipate  that  our  semlclasslcal  method  will  be  extremely  useful  for 
this  purpose. 
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DISCUSSION  OF  THE  EFFECTS  OF  CLUSTERING  ON  ION  STABILITY 
AND  NEUTRALIZATION  RATES 


We  consider  our  discovery  of  a  closed-form  perturbation  solution 
for  the  general  scattering  problem  for  the  lon-dlpole  Ion  system, 
Including  capture  collisions,  to  be  of  the  greatest  Importance.  Not 
only  does  It  allow  rapid  and  economical  estimation  of  capture  cross 
section  and  rates,  but  even  more  Important,  it  makes  possible  a 
methodical  investigation  of  the  effect  of  changing  conditions  and 
parameters.  As  our  experience  with  the  properties  of  this  solved  problem 
grows  we  expect  our  understanding  to  become  much  more  complete  and 
secure,  and  we  will  then  be  able  to  extend  It  to  the  related  but 
somewhat  more  complex  problems  that  arise  in  considering  the  neutraliza¬ 
tion  of  clustered  Ions. 

The  commonest  clustering  species  associated  with  most  ions  in  the 
lower  Ionosphere,  both  positive  and  negative,  is  water.  The  most 
Important  clustered  Ions,  therefore,  have  highly  polar  components 
(water  molecules)  very  much  subject  to  the  effects  of  transient 
electrostatic  forces  caused  by  a  passing  Ion  of  opposite  sign.  In  most 
cases  the  clustered  Ions  are  expected  to  be  very  soft,  with  bending 
modes  of  low  frequency  exceptionally  easily  excited,  since  the  lon- 
dlpole  attraction  that  creates  the  clustering  Is  generally  less  strongly 
directionally  dependent  than  most  chemical  bonds.  Furthermore,  a 
comparatively  simple  electrostatic  model  for  the  structure  allows 
reasonable.  If  rough,  estimation  of  the  restoring  forces  and  the 
characteristic  bending  frequencies.  Information  of  this  type  Is  needed 
to  develop  a  more  quantitative  approach  to  the  neutralization  rates  of 
these  clustered  species. 
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As  part  of  our  progreun  we  Intend  to  extend  our  approach  from 
rotational  excitation  of  diatomic  molecules  to  the  excitation  of 
bending  modes,  l.e..  Internal  rotations  In  a  simple  hydrated  Ion 
(for  Instance,  an  Ionic  core  such  as  an  alkali  Ion  with  a  single 
attached  water  molecule). 

Generally,  hydration  will  have  three  effects  on  an  Isolated 
molecular  Ion:  It  will  Increase  the  stability  of  the  Ion  (whether 
It  Is  positive  or  negative);  It  will  Increase  the  number  of  degrees 
of  freedom  of  the  molecule;  and  especially  It  will  produce  low 
frequency  Internal  modes  (soft  bending  vibrations).  Increasing  the 
stability  of  the  negative  Ion  contracts  the  exponential  tall  of  the 
electron  cloud  and  causes  more  rapid  fall-off  with  distance  of  the 
overlap  matrix  leading  generally  to  a  lower  neutralization 

cross  section.  At  the  same  time.  Increased  stability  of  the  Ions 
relative  to  the  neutrals  causes  all  the  neutral  molecular  potential 
curves  of  the  combined  molecule  to  move  upward  relative  to  the 
coulomblc  Ion- Ion  potentials,  and  causes  each  crossing  point  to  move 
outward  to  larger  R.  In  general,  therefore,  hydration  tends  to 
reduce  the  cross  section  for  neutralization  by  electron  capture 
(as  described  by  the  curve -crossing  model).  On  the  other  hand,  by 
Increasing  the  number  of  Internal  degrees  of  freedom,  and  especially 
of  the  soft  bending  modes  that  can  easily  be  excited  by  the  tidal 
forces  associated  with  the  close  passage  of  a  second  Ion  of  opposite 
charge.  Increasing  hydration  Is  associated  with  an  Increasing  cross 
section  for  neutralization  by  tidal  capture.  This  Is  a  to^ally 
different  mechanism  from  electron  transfer  (although  temporary  capture 
can  sometimes  also  allow  multiple  passages  past  the  crucial  curve- 
crossings,  thus  enhancing  the  opportunities  for  electron  transfer). 
Thus,  In  general,  we  can  envisage  the  possibility  that  hydration  may 
have  two  competing  effects  on  Ionic  neutralization  rates:  It  decreases 
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the  electron-transfer  rate,  and  it  ultimately  compensates  by  increasing 
the  tidal  recombination  rate.  In  some  cases,  and  perhaps  in  many,  the 
combined  effect  may  lead  to  a  minimum  in  the  recombination  rate  as  a 
function  of  n,  the  degree  of  hydration  of  one  of  the  ions,  at  some  small 
value,  perhaps  n  >  1  or  2. 

These  effects  can  be  enhanced  further  by  the  occurrence  of  special 
structural  configurations.  We  have,  therefore,  been  led  to  consider 
the  structures  associated  with  hydration  in  positive  and  negative  ions. 

In  the  case  of  water  attached  to  a  positive  ion,  the  principal 

binding  can  be  considered  to  arise  from  simple  ion-dipole  binding. 

In  the  case  of  hydrated  negative  ions,  on  the  other  hand,  it  seems 

likely  that  the  attachment  of  a  water  molecule  will  be  centered  at  one 

of  the  H  atoms,  and  that  a  form  of  hydrogen  bond  is  likely  to  result. 

Ordinary  hydrogen  bonds  tend  to  be  comparatively  straight,  with  the 

H  atom  roughly  on  the  line  of  centers  between  its  two  neighbors,  but 

usually  asymmetrically  placed  much  closer  to  the  atom  with  which  it  has 

its  primary  chemical  bonding.  Under  very  exceptional  conditions, 

however,  symmetrical  hydrogen  bonds  have  recently  been  observed. 

Pauling^^  estimates  the  typical  restoring  constant  for  the  bending  of 

2 

a  hydrogen  bond  in  water  to  be  0.003  kcal/mole-deg  .  Whether  this  is 
a  typical  figure  that  can  be  used  for  estimating  bending  force  constants 
in  hydrated  negative  ions  is  uncertain,  but  it  at  least  provides  some 
guidance . 

In  some  cases,  exceptional  structures  may  exist  that  will  provide 
special  stability  to  certain  ions,  and  cause  various  parameters, 
including  bending  force  constants,  to  deviate  from  the  general  pattern. 
As  an  example,  let  us  consider  the  special  case  of  the  ion  NO2  'H^O. 

This  is  probably  the  principal  ion  of  mass  64  observed  in  a  decaying 
plasma  in  air  by  Dr.  Merle  Hirsch  (reported  at  the  28th  Gaseous 
Electronics  Conference  on  October  21-24,  1975).  This  ion  apparently 
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recombines  exceedingly  slowly,  with  a  rate  coefficient  of  about 
"8  2 

10  cm  /sec,  and  appears  to  have  unusually  high  mobility  for  a 
negative  ion.  It  would  be  of  considerable  Interest  to  make  further 
measurements  on  this  Ion  and  even  to  determine  Its  structure.  Several 
possible  structures  can  be  considered. 

Structure  A,  Illustrated  In  Fig.  2,  appears  to  be  a  plausible 
chemically  bonded  form  for  what  may  be  the  lowest  stable  state  of  the 
ion  that  case,  the  three  oxygen  atoms  may  be  bonded  to 

nitrogen  roughly  as  in  NK^  (i.e.,  in  a  pyramidal  structure).  This 
structure  is  not  a  simple  hydrate,  but  requires  considerable  rearrange¬ 
ment,  and  if  it  is  produced  at  all,  it  is,  no  doubt,  the  result  of  a 
more  complicated  series  of  steps  than  simple  hydration.  It  might, 
therefore,  be  more  common  under  high-pressure  conditions  than  at  low. 

Structure  B  represents  a  typical  hydrate  of  NO^  ,  assuming  a 
standard,  almost  linear,  hydrogen  bond.  Such  a  structure  is  rather 
open  and  loose  and  subject  to  fairly  free  rotation  of  the  components 
about  the  H-bond,  but  electrostatic  forces  arising  from  the  contributing 
structure  B'  may  make  the  coplanar  structure  somewhat  more  stable. 

One  can  also  speculate  on  the  possibility  of  a  more  symmetric 
coplanar  structure  as  illustrated  in  C,  which  might  also  have 
significant  contributions  of  the  companion  structure  D.  These 
structures  have  the  possibility  of  additional  stabilization  through 
delocalization  of  the  electron  and  formation  of  a  second  hydrogen  bond, 
but  this  is  achieved  at  the  expense  of  bending  the  hydrogen  bond  and 
possibly  by  forcing  the  hydrogen  to  a  more  symmetric  position.  If  the 
structure  C  (and  D)  exists,  it  is  sure  to  be  coplanar  and  rather  stiff 
in  its  bending  modes,  and  this  will  have  consequences  for  both  the 
energetics  and  the  entropy  of  the  molecule.  A  theoretical  examination 
of  these  possible  structures  as  well  as  further  experimental  studies 
of  the  properties  of  this  hydrated  ion  would  be  well  warranted.  If 
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FIGURE  2  STRUCTURES:  NOj’.HjO  =  H^NOj'  (MASS=64) 
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it  has  a  compact  structure  like  either  A  or  C,  its  exceptionally  low 
recombination  rate  might  be  understood. 

Once  the  ion  is  doubly  hydrated,  a  much  looser  structure  generally 
results,  with  very  significant  softening  of  the  bending  modes  as  well 
as  other  consequences  such  as  increased  size  and  lower  mobility.  Among 
the  possibilities  are  such  structures  as  E  and  F,  of  which  E  is 
obviously  derived  from  A  and  F  from  B  or  C. 

If  the  ring  structures  we  have  discussed  here  are  valid  in 
NO2  *1120,  they  may  also  apply  to  other  negative  ions  such  as 
NO^  and  CO^  when  they  are  singly  hydrated.  In  that  case  a  number  of 
important  singly  hydrated  ioiiS  may  be  quite  different  in  their 
properties  from  the  higher  hydrates.  These  effects  would  appear  both 
in  the  thermodynamics  of  the  ions  and  their  spectroscopy,  and  also  in 
other  properties  such  as  neutralization  rates,  photodissoci  tion  and 
photodetachment,  and  possibly  in  such  collision  properties  as  mobility. 

These  arguments  suggest  that  singly  hydrated  ions  will  often  have 
rather  special  properties.  Their  electrons  will  be  more  strongly 
bound  than  those  of  the  parent  unhydrated  ions,  which  will  tend  to 
result  in  lower  neutralization  rates  by  the  electron  transfer  (curve¬ 
crossing)  mechanism.  On  the  other  hand,  having  stiff  and  compact 
structures,  they  will  also  have  rather  low  cross  sections  and  rates 
for  tidal  capture  into  an  elliptical  orbit  by  ions  of  opposite  charge, 
the  process  that  we  expect  to  be  rate-determining  for  recombination  of 
large  cluster-ions  at  low  pressure. 
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A  uoiforir  Bemiclaaaical  S  matrix  haa  been  developed  for  colliaiona  of  charged  parti- 
clea  with  rolttiog  rigid  dipolea,  with  uae  of  first-order  perturbation  theory.  The  re¬ 
sulting  expression  is  analytical,  depending  on  tidhulated  functions,  and  trivial  to  calcu¬ 
late;  it  allows  evaluations  of  quantum  transitions  in  classically  forbidden  regions,  and 
of  quantum  interference  effects. 


In  the  course  of  studies  of  the  scattering  of 
electrons  and  ions  by  simple  polar  molecular 
targets,  we  have  discovered  a  versatile  analyti¬ 
cal  form  for  the  S  matrix  in  the  limit  of  semi¬ 
classical  perturbation  theory  when  the  anisotrop¬ 
ic  part  of  the  interaction  is  dominated  by  a  di¬ 
pole  term. 

Cross'  developed  a  classical  perturbation  theo¬ 
ry  of  ion-molecule,  including  ion-dipole,  scatter¬ 
ing  in  a  formulation  emphasizing  angular  coordi¬ 
nates  and  spherical  trigonometry.  We  find  it  ad¬ 
vantageous  to  make  maximum  use  of  angular  mo¬ 
mentum  conservation,  formulating  the  problem 
in  the  classical  version’*  ’  of  a  (J,  M,j,  1)  angular 
momentum  coupling,  where  j  is  the  angular  mo¬ 
mentum  of  the  target  molecule  and  /  is  the  col- 
lisional  angular  momentum,  while  J  and  .Vf  refer 
to  the  total  angular  momentum  and  its  projection. 
[Where  the  distinction  is  important,  j,  I,  and  J 
will  be  used  for  quantum  numbers,  and  the  cor¬ 
responding  classical  angular  momenta  are  given 
by  expressions  like  j''  =  {j  +  ^)/).|  We  have  used 
the  semiclassical  5-matrix  formulation  of  Mil¬ 
ler,’  combined  with  first-order  perturbation  dy¬ 
namics.  By  applying  a  canonical  Hamiltonian 
transformation,  we  find  the  phase  of  the  5  matrix 
to  have  a  contribution  not  identified  by  Cross. 

The  uniform  semiclassical  form  for  the  5  ma¬ 
trix  provides  an  appropriate  estimate  of  quantum 
effects  in  the  nonclassical  tunneling  region  and 
at  the  classical  turning  points,  as  well  as  inter¬ 
ference  effects  in  the  classically  allowed  region. 
In  addition,  by  use  of  the  semiclassical  5  matrix 
in  (J,  M,j,  1)  coupling  followed  by  quantal  recou¬ 
pling  using  6-j  symbols  and  construction  of  the 
scattering  amplitude  using  the  Wigner  rotation 
matrices,  we  can  obtain  the  quantal  diffraction 
effects  in  small-angle  scattering. 

We  consider  a  problem  with  a  zero-order  po¬ 


tential  interaction  in  which  j* and  f  are  separately 
conserved: 

Vo=Vo(R)-  (1) 


R  is  the  radial  coordinate  associated  with  the 
collision.  In  the  plane  determined  by  T,  the  un¬ 
perturbed  collisional  motion  follows  a  classical 
trajectory  given  by  R(t)  and  ♦(/),  and  in  the  plane 
determined  by  |  the  rotor’s  unperturbed  motion 
is  given  by  e(/).  The  angle  k  between  the  two 
planes  is  fixed  by  the  magnitudes  j,  I,  and  J: 


N-cosk  = 


I*T 

Tiitr 


(2) 


The  line  of  intersection  of  the  two  planes — or  the 
plane  perpendicular  to  it — defines  a  reference 
direction  for  the  measurement  of  4>  and  6,  each 
in  its  own  plane.  It  is  convenient  to  take  /  =  0  at 
a  central  point  in  the  trajectory,  and  assume  the 
unperturbed  R{t)  to  be  an  even  function  of  t.  We 
can  then  write  0(/)  =  0o  +  S(f)  and  =  + 

where  Q(l)  and  $(/)  are  odd  functions  of  t. 

The  anisotropic  interaction  will  be  treated  as 
a  perturbation  (to  which  an  isotropic  part  can  be 
added  as  well).  Particular  simplifications  occur 
for  a  dipole  interaction,  and  our  attention  will  be 
limited  to  that  case;  extensions  to  higher  multi¬ 
poles  and  nonrigid  targets  are  possible,  but  they 
will  probably  entail  greater  relia.'ce  on  numeri¬ 
cal  methods.  We  shall  assume  the  general  form 
of  the  interaction  to  be 


V,{R,t)  =  U{R)cosr,  (3) 

where 


COS7  =  -  cos*  COS0  +A'sin*  sine 
=  - 1(1  +N)cos4jl) 

-^(l-iV)cosi/.(/),  (4) 
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with 

0.(/)=e(/)±4>(/)=(i.''+$,(/),  (5) 

and 

i/./=0o^*o.  fjn=eu)t$(/).  (6) 

The  dynamic  effects  of  the  perturbation  can 
now  be  evaluated.  The  associated  action  integral 
is  a  function  of  Og  and 

A  =  A(eo,<ho)  =  COS} (l))dt.  (7) 

This  is  conveniently  simplified  through  use  of  the 
angular  variables  and 


A  =A(i^ /,*.“) 

=  -i(l  +N)B*  cosi/," 

-|(1-N)B.  oosi^.“. 

(8) 

where 

B,=JlU(RU))  cos^,(/)<i'(. 

(9) 

Solving  the  Hamiltonian  equations  of  motion  for 
the  changes  in  /  and  J  to  first  order,  one  can 
easily  show  that 


300 

If  we  define 

A*,  j  lA/), 
and 

z,  =  (2fir‘(uy)B,. 
we  have  the  simple  result 

A*,  =  -  ;r  ‘  dA/dil  -  z .  sin^l  ° . 


(10) 

(11) 

(12) 

(13) 


The  range  of  classically  accessible  values  of  Ak^ 
and  Ak.  is  determined  by  Z,,  which  depends  on 
the  unperturbed  values  of  /  and  i,  on  J,  and  on 
the  collisional  energy  K,  In  addition,  each  value 
of  Ak^  or  Ak.  is  associated  with  a  pair  of  per¬ 
missible  values  for  the  associated  angle  4^°  or 
4.°,  so  each  set  of  values  {Ak\,  Ak. )  can  be 
reached  by  four  different  classical  trajectories. 
These  will  lead  to  an  interference  pattern  in  the 
S  matrix. 

The  classical  phase  is  expressed  by 


ip  =  -fdl[R  '(/)p'(/)  +♦'(/)/'''(/)  +0 '«)  J'"‘(/)l 


to  the  original  variables  R,  p,  etc.,  through  a 
canonical  transformation)  chosen  in  such  a  way 
that  the  coordinates  R  ,  Q',  and  appear  ex¬ 
plicitly  only  in  the  perturbation  potential  term 
(Y,)  in  the  transformed  Hamiltonian.  With  the 
help  of  the  appropriate  canonical  transformation, 
it  can  be  shown  that 

,p  =  -Q^Aj^*JiAl-A(Qo,  *o) 

^-4AAk,-^.°hAk.  -AU.\4.°).  (15) 

Making  the  appropriate  substituti')ns,  we  can 
write  this  in  the  form 

(P/fi  -t-T}.  ,  (16) 

where 


T}^=-tl^°Ak^+Z,COS4° 

=  (Z/- Ait ,*)■''* +  sin’ ’(Air, /Z.).  (17) 


To  construct  the  S  matrix,  we  need  the  Jacob¬ 
ian  D  as  well  as  the  phase  ip: 

8(7,/)  d{Aj,Al) 

~  8(00,  4.0)  ■  8(00,  *0) 


if  we  take 


,  dAk\  ^Ak. 

dc.o 


=c,c. , 


(18) 


C.=-2^^  =  2Z,cosC.“ 

8''  . 

=  2(Z.'-A<?. *)■'*.  (19) 

The  symmetry  of  the  solutions  is  such  that  noth¬ 
ing  is  changed  if  we  replace  (Gq,  ^o^  by  (0o  +  Jr, 
♦0*  ’f),  which  leads  to  a  selection  rule  of  Ak\ 

=  integer. 

The  ii  matrix  is  the  sum  of  four  terms,  corre¬ 
sponding  to  the  four  pairs  of  roots  (. .“,  ».°)  al¬ 
lowed  by  Eq.  (13),  each  being  of  the  form 

[(2;7/)'/)1  ■ '  *exp(jV/f/) .  (20) 

Exercising  appropriate  care  to  identify  the  phas¬ 
es,  and  making  use  of  the  form  of  D  and  v’,  we 
find  that  we  can  combine  terms  so  that  S  can  be 
factored  into  two  independent  portions, 

S=^,.S.,  (21) 

where 

i,  =(4/7rC,)''*sin(U+(3,)  (22) 


+  4>'(/) 


BV, 

8*' 


(14) 


where  R',  p',  etc.,  are  “new"  variables  (related 


with 

+ cos'‘(-Ai8/Z,).  (23) 

This  form  is  appropriate  well  into  the  classical 
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region  of  motion,  orC/»0.  Quan¬ 

tum  effects  become  important  near  C,  =0,  and 
there  we  can  make  use  of  a  uniform  approxima¬ 
tion  in  the  Bessel  function  form,^’  ^  which  re¬ 
duces  to  (22)  in  the  limit  of  large  C,: 

(24) 

When  C,  becomes  imaginary,  (24)  goes  over  into 
an  exponential  tunneling  form  valid  for  Ait/ 

» Z/,  or  C/« 0: 

(25) 

The  form  (24),  however,  is  valid  everywhere, 
and  is  essential  near  the  edge  of  the  classical 
region,  where  C/-0,  i,e,,  at  the  classical 
turning  points. 

The  S  matrix  is  desired  in  terms  of  the  initial 
and  final  angular  momenta,  which  are  connected 
with  .;,  /,  A  j,  and  A/ by 

/,  =  /-A//2,  /j  =  /  +  A//2, 

/,=7-a//2,  j2=.;-^Aj/2. 

Also  the  initial  and  final  collisional  energies  are 
connected  with  the  average  value  E  by 

E*{j  +  h  )-H "  /2/  =  £  1  +  ( ,/ ,  +  i  )Vi  *  /2I 

=  £,  +  (./,  +  ^)-//V2/.  (27) 

With  use  of  these  relations  [(26)  and  (27)J  and  the 
properties  of  the  Bessel  functions,  it  is  easy  to 
show  that  this  ii  matrix  is  correctly  symmetric 
and  approximately  unitary — approaching  unitari- 
ty  correctly  as  I  and  j  grow  large.  It  is  not  diffi¬ 
cult  to  renormalize  it  for  small  j  and  I  to  make 
it  correctly  unitary. 

We  have  evaluated  the  integrals  5,  [Eq.  (9)|  for 
two  specific  cases,  in  both  of  which  the  perturba¬ 
tion  is  an  ion-  or  electron-dipole  interaction, 

r,  =(fic 'R^)c0S7  .  (28) 

When  the  target  is  a  neutral  dipole,  case  (a),  the 
unperturbed  potential  vanishes,  To  =  0.  The  re¬ 
sult  can  be  expressed  in  terms  of  a  parameter 
measuring  the  ratio  of  angular  velocities  of  the 
rotor  and  the  collisional  motion  at  /  =0, 

=(j  +  2)(/  +  i)//- '2/£,  (29) 

where  1  is  the  moment  of  inertia  of  the  dipole,  a 
factor 

nuic'tiHl-l),  (30) 

where  w  is  the  reduced  mass  in  the  collision. 


and  the  integrals 

=a.,LA:,(u)/^/foK)J.  (31) 

The  same  integrals  appear  in  the  problem  as  for¬ 
mulated  by  Cross,'  and  also  in  Percival’s  theory 
of  excitation  of  hydrogenic  atoms  by  charged  par¬ 
ticles.®  Then  the  quantity  Z,  becomes 

Z,  =  (wMe///*)(/  +  i)-(1±Ar)/,.,(wJ.  (32) 

We  are  using  this  solution  in  a  study  of  scattering 
and  rotational  excitation  in  collisions  of  electrons 
or  ions  with  polar  molecules. 

When  the  dipolar  target  has  a  net  charge,  and 
the  bombarding  species  is  an  ion,  is  the  Cou¬ 
lomb  potential,  Va  =  ±e^/R.  We  have  found  an 
interesting  solution  for  a  special  case,  where  the 
unperturbed  collisional  energy  £  is  0  and  the  tra¬ 
jectory  is  parabolic.  In  that  case,  (b),  we  have 

Wb  =  (.;  +  |)a  +  i)'«VwcV,  (33) 

and  the  integrals  7,  ,,  lead  to  Airy  functions, 

^i.b  -27rL+.rAi(.r)-Ar*^*Ai'(x)j, 

x  =  (34) 

This  solution  is  of  interest  in  connection  with  cap¬ 
ture  or  detachment  between  free  (hyperbolic)  tra¬ 
jectories  with  eccentricity  slightly  greater  than 
1,  and  large  elliptic  orbits  with  eccentricity 
slightly  less  than  1. 

The  separation  of  the  equations  connecting  the 
phase  angles  i/  /  with  the  momenta  Afe,  in  the 
simple  form  of  Eq.  (13)  appears  to  be  peculiar  to 
the  dipole  angular  dependence  of  K,.  It  is  inde¬ 
pendent  of  the  R  dependence  of  and  and  so 
the  factorization  of  the  5  matrix  and  the  general 
form  (24)  of  its  factors  persists,  only  the  inte¬ 
grals  or  7 ,  changing  with  the  functional  form 
of  the  7?  dependence.  These  integrals  can  be  ob¬ 
tained  by  numerical  quadrature  when  an  analyti¬ 
cal  form  is  not  available. 

When  the  angular  dependence  of  V,  is  not  that 
of  a  simple  dipole,  P^(,cos■^),  the  elimination  of 
the  phase  angles  and  the  expression  of  S  as  a 
function  of  the  angular  momenta  and  their  per¬ 
turbations  is  not  so  simple  as  in  the  dipole  case, 
and  may  have  to  be  carried  out  by  numerical  in¬ 
terpolation.  In  that  case,  the  S  matrix  will  prob¬ 
ably  not  factorize. 
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APPENDIX  B 

Electron  loss  cross  sections  for  0~,  O2 »  NO^,  and  NO^  in 
several  gases* 

R  A.  Bennett^  J.  T.  Moseley,  and  J.  R.  Peterson 

Stanford  Research  Institute,  Menlo  Park,  California  94025 
(Received  18  November  1974) 

The  croM  sections  for  electron  lou  for  0‘.  Of.  NOf,  and  NOf  on  He,  Nj,  and  Ar,  for  O"  on  O2. 
and  NOf  on  NO]  have  been  measured  using  a  beam  attenuation  technique.  The  measured  cross 
sections  are  the  sum  of  coUisional  detach.nent  and  electron  transfer  cross  sections.  Except  for  O'  + 

O],  the  cross  sections  are  nearly  constant  for  ion  energies  in  the  range  1  to  4  keV  and  generally 
increase  with  target  mast. 


INTRODUCTION 

We  report  the  measurement  of  electron  loss  cross 
sections  lor  O',  Oi,  NOi,  and  NOj  In  several  gases. 

The  work  originate  In  an  attempt  to  determine  any  ex¬ 
cited  state  fractions  of  negative  Ion  beams  used  In  this 
laboratory  for  lon-lon  recombination  rate  measure¬ 
ments.*  A  beam  attenuation  technique  similar  to  that 
of  Turner  et  was  used.  No  evidence  for  suspected 
excited  state  ions  was  found,  but  the  e}q>erimental  setup 
made  It  possible  to  determine  electron  loss  cross  sec¬ 
tions.  The  total  loss  cross  section  which  is  here  asso¬ 
ciated  primarily  with  colllslonal  detachment,  may  be 
determined  by  measuring  the  current  /  remaining  when 
an  Ion  beam  of  Initial  current  /g  traverses  a  distance  I 
through  a  chamber  containing  a  gas  with  number  density 
n.  If  the  beam  contains  two  components  (such  as  a 
ground  state  and  an  excited  state)  in  fractional  amounts 
/i  and  /{ =  1  - /i,  with  different  loss  cross  sections  Q|  and 
Qt ,  respectively,  the  current  collected  will  be 

+  •  (1) 

A  plot  of  log  (///g)  versus  pressure  will  thus  show  a 
curve  with  different  low-pressure  and  high-pressure 
asymptotic  slopes.  The  two  slopes  determine  the  two 
characteristic  cross  sections  Qj  and  Q,,  and  the  zero 
pressure  intercepts  yield  /,  and  /g.  For  every  case  re¬ 
ported  here,  the  attentuatlon  curve  of  log  I/Iq  versus 
pressure  appears  to  be  a  single  straight  line,  indicating 
a  single  state  for  the  ions  in  the  beam,  or  at  least  that 
the  products  fQ  lor  all  components  were  the  same  to 
within  about  5%.  The  possiblility  of  loss  mechanisms 
other  than  colllslonal  detachment  will  be  discussed  below 

APPARATUS 

The  apparatus  Is  a  merged  beams  system*  which  was 
modified  to  make  attenuation  measurements  as  shown  in 
Fig.  1.  Negative  ions  extracted  from  a  duoplasmatron 
ion  source  were  accelerated  to  the  beam  energy  of  1  to 
4  keV.  After  mass  selection  by  a  crossed  fields  velocity 
filter,  the  ions  were  focused  and  deflected  through  a 
1  cm  aperature  into  an  attenuation  chamber.  The  fo¬ 
cused  beams  were  usually  less  than  0. 6  cm  In  diameter, 
so  that  losses  to  the  aperture  boundaries  were  negligible. 
After  traveling  a  length  f  -  34. 3  cm,  the  beam  was  col¬ 
lected  in  a  shielded  Faraday  cup  and  monitored  by  an 
electrometer.  The  secondary  electron  suppression  grid 
In  front  of  the  Faraday  cup  was  maintained  at  -  30  V, 


thus  also  preventing  collection  of  slow  electrons  that 
were  colllslonally  detached  in  the  target  gas  along  the 
beam  path.  Tests  using  a  two-part  concentric  Faraday 
cup  established  that  scattering  losses  of  the  beam  were 
insignificant.  Differential  pumping  limited  the  measured 
attenuation  to  the  last  chamber. 

The  attenuation  chamber  pressure  was  measured  using 
a  Plranl  gauge  and  a  Varlan  nude  Ion  gauge,  each  of 
which  was  calibrated  against  an  MKS  Baratron  capaci¬ 
tance  manometer  down  to  10“^  Torr.  Both  gauges  were 
observed  to  be  quite  linear  with  pressure  above  10*^ 
Torr,  and  were  assumed  to  be  linear  down  to  10"*  Torr, 
The  ion  gauge  was  used  for  the  pressure  measurements 
on  the  rare  gases  and  onNg,  and  the  Plranl  gauge  was 
used  for  Og  and  NOg ,  because  for  these  gases  stability 
problems  were  encountered  when  using  the  ion  gauge. 
However,  some  attenuation  measurements  were  made 
on  each  gas  using  both  gauges. 

As  a  test  of  the  experimental  technique,  the  attenua¬ 
tion  by  argon  of  3  keV  NO*  ions  from  a  duoplasmatron 
source  was  studied.  The  attenuation  curve  showed  a 
two-slope  behavior  similar  to  that  seen  by  Mathis  et 
ah*  Most  of  the  NO*  beam  had  an  attenuation  cross  sec¬ 
tion  of  16  X  10'**  cm*.  A  small  fraction  of  the  NO*  beam, 
about  15%,  had  a  much  larger  cross  section.  The  exact 
fraction  depended  on  the  ion  source  parameters,  includ¬ 
ing  the  gas  mixture  and  the  pressure. 

With  the  use  of  various  mixtures  of  Og,  NO,  and  NOg 
in  the  negative  ion  source,  beams  of  10**  A  of  O’,  10** 

A  of  Og,  10’*  A  of  NOg,  and  lO"*  A  of  KOj  could  be  ob¬ 
tained.  The  ion  source  pressure  was  in  the  range  of 
10  to  100  mTorr.  Only  long-lived  ion  states  contributed 
to  the  attenuation  measurements,  because  the  minimum 
flight  time  from  ion  source  to  attenuation  ci>amber  is 
7  fisec  with  4  keV  O'  ions. 

RESULTS 

The  attenuation  of  a  3. 75  keV  beam  of  Og  ions  by  ni¬ 
trogen  is  shown  in  Fig.  2.  The  logarithm  of  the  col¬ 
lected  current  was  plotted  versus  the  nitrogen  pressure 
indicated  by  an  Ion  gauge.  The  residuals  of  a  straight 
line,  least  squares  fitted  to  the  data,  showed  no  system¬ 
atic  variation.  For  the  same  beam  energy  and  target, 
the  linearity  and  slope  of  the  attenuation  curve  were  very 
reproducible.  The  slope  of  the  attenuation  curve  as 
plotted  in  Fig.  2  yielded  a  cross  section  of  12. 8  ±  1. 2 
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FIG.  1,  Exparimeatal  appanitus. 


X  10*‘*  cm*  tor  the  attenuation  of  3. 75  keV  OJ  by  nitrogen, 
using  Eq.  (1)  with  /t  - 1  and  after  corrections  to  the  gauge 
pressure  readings.  The  uncertainty  In  the  absolute  mag¬ 
nitude  of  the  cross  section  so  determined  arose  from  the 
uncertainties  In  current  and  pressure  measurement.  The 
uncertainty  In  the  current  measurement  was  taken  to  be 
twice  the  rated  electrometer  accuracy,  or  4%.  The  un¬ 
certainty  In  the  pressure  measurement  was  conservative¬ 
ly  placed  at  ±15%  In  He  and  0,,  ±13%  In  Ar,  and  ±5% 


N,  PRESSURE  (1(r*T) 
FIG.  2.  Attenaatton  of  0{  by  N]. 


In  Nt  and  NOj. 

The  attenuation  cross  sections  for  O',  Of,  NOf,  and 
NOf  In  various  gases  as  functions  of  laboratory  ion  ener¬ 
gy  are  shown  In  Figs.  3-6.  The  cross  sections,  at  a 
representative  3  keV  Ion  energy,  are  tabulated  In  Table  I. 
In  He,  0| ,  and  NO^ ,  which  may  form  negative  Ions,  the 
electron  loss  cross  section  Is  the  sum  of  cross  sections 
for  charge  exchange  and  colllslonal  detachment. 

In  addition  to  single  electron  loss  collisions,  with 
whose  cross  sections  a.,g  we  are  primarily  concerned 
here,  we  must  consider  two-electron  loss  collisions 
which  form  positive  Ions  In  the  beam  and  subtract  from 
the  negative  Ion  current.  Thus,  as  Rlsley  and  Geballe* 
have  pointed  out,  these  negative  Ion  attenuation  measure- 


LAB  ENERGY  (keV) 
FIG.  3.  Electron  loss  cross  sections  for  O'. 
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TABLE  I,  Electron  lose  cross  sections 
for  O',  C^,  no;,  and  NO;  at  3  keV. 


O' 

o<10‘’*cm*) 

oi  no; 

NCS 

He 

4.6 

4.4 

s.o 

7.0 

Ng 

10.7 

13.1 

16.0 

15.5 

Oj 

24 

ess 

ess 

s  s  • 

Ar 

12.5 

21.8 

16.0 

18.3 

NO, 

•  S  • 

ess 

•  •  s 

39.5 

ments  In  principle  measure  single  collision  cross  sections 
containing  u.io+Zo.n,  and  at  high  energies  one  should 
consider  even  higher  ionization  processes.  In  the  present 
case,  however,  o.„  is  probably  at  least  20  times  smaller 
than  (T.jg ,  judging  trom  the  low  energy  trends  in  the 
cross  sections  obtained  by  Matid  and  Cobid.’  We  there¬ 
fore  consider  these  measurements  to  represent  primar¬ 
ily  single  electron  loss  cross  sections  o.jg .  In  most 
cases  the  loss  mechanism  is  colllslonal  detachment. 
However,  in  the  cases  of  Og  and  NOg  targets  which  have 
positive  electron  affinities,  charge  transfer  may  also 
be  an  Important  process.  Unfortunately  only  one  mea¬ 
surement  was  made  in  each  of  these  gases.  Charge 
transfer  to  He,  which  has  an  excited  negative  ion  state, 
is  about  19  eV  endothermic  and  is  almost  certainly  neg¬ 
ligible.  Similarly,  excitation  of  the  2  eV  excited  virtual 
state  in  Ng'  is  probably  unimportant  as  a  loss  chaimel. 
Some  of  these  effects  will  be  considered  below. 

DISCUSSION 

In  Fig.  7  the  attenuation  cross  sections  for  O*  in  He 
and  Ar  are  plotted  against  O*  velocity  and  can  be  com¬ 
pared  with  other  results.  The  He  results  are  in  quite 
good  agreement  with  Hasted’s  results'  which  also  agree 
with  the  low  energy  data  given  by  Wynn,  Martin,  and 
Bailey  (WMB).’  No  uncertainty  estimates  are  given  by 
Hasted,  and  it  is  not  ^parent  that  the  15%  difference  in 
results  is  significant.  The  relative  consistency  among 
our  data  points  indicates  that  there  is  little,  if  any,  en¬ 
ergy  dependence  in  the  cross  section  at  these  energies. 

The  O*  +  Ar  results  of  Hasted  are  about  30%  lower  than 
ours  and  again  match  the  WMB  data  at  low  energy.  The  re¬ 
sult  B  of  Roche  and  Goodyear*  are  slightly  below  those  of 


LAB  ENERGY  (keV) 

FIG.  5.  Electron  loss  cross  sections  for  NOj. 


WMB.  Between  the  speeds  of  lOand  22x  10*  cm/sec,  the 
Hasted  results  for  both  He  and  Ar  show  a  slight  increase, 
whlleourcrosssectionsareessentlallyconstant.  The  re¬ 
sults  of  Matic  and  Coble*  on  O' -t-  Ar  begin  at  8  keV,  and 
therefore  it  is  not  possible  to  compare  them  directly 
with  our  results  or  those  of  Hasted.  However,  if  the 
results  at  4  and  8  keV  are  both  correct  in  absolute  mag¬ 
nitude,  the  cross  section  must  exhibit  an  unexpectedly 
strong  energy  dependence  in  the  energy  region  between 
the  measurements. 

A  comparison  with  other  results  for  O'  in  Ng  and  Og  is 
made  in  Fig.  8.  The  results  of  Hasted  and  Smith*  for 
O'-^Ng  are  about  20%  below  ours  and  show  a  sL'ght  de¬ 
crease  with  speed  while  our  results  show  a  slight  in¬ 
crease.  Apparently  the  O*  Ng  cross  section  is  about 
constant  in  this  energy  range.  The  results  of  Matic  and 
^obi£  for  O*  *  Ng  extend  down  to  5  keV  and  again  are 
substantially  below  our  results  at  4  keV. 

The  O'-f  Og  data  are  in  reasonably  good  agreement 
with  Hasted  and  Smith  in  both  magnitude  and  energy  de¬ 
pendence.  The  results  of  Bailey  and  Mahadevan**  and 
Ranjan  and  Goodyear**  are  somewhat  lower.  The  cross 
section  is  larger  than  would  be  expected  from  the  relative 
size  of  Og  compared  to  Ng  and  Ar.  The  energy  depen¬ 
dence  in  the  1  -4  ke  V  range  is  in  distinct  contrast  to  that 
for  He,  Ar,  and  Ng.  This  can  be  explained  by  the  open 
charge  transfer  channel,  O'-fOg-O-fOg,  which  supple¬ 
ments  colllslonal  detachment  as  a  loss  mechanism.  The 
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FIG.  7.  O'  electron  loss  cross  sections  In  He  and  Ar.  TW 
=  this  work.  H  =  Hasted.''  WMB  =  Wynn,  Martin,  and  Bailey. 
RG  =  Hocheand  Goodyear,*  MC  =  Mati6  and  Cobid.* 


charge  transfer  data  of  Mathis  and  Snow**  for  O'  +  O*  are 
similar  to  those  of  Snow,  Rundel,  and  Geballe**  and  are 
indicated  in  Fig.  8.  These  data  suggest  that  the  colli - 
sional  detachment  cross  section  for  O'-t-Oj  is  about  con¬ 
stant  it  15  Y 10'**  cm*  between  1  and  4  keV. 

The  only  reported  electron  loss  measurements  for 
Oj  are  those  by  WMB’  for  the  case  Oi+  He.  Their 
results  reach  only  as  high  as  400  eV,  but  the  ^^iparent 
asymptotic  limit  of  their  data  agrees  well  with  what  we 
find  at  higher  energies. 

We  know  of  no  previous  measurements  of  electron 
loss  cross  sections  for  NOi  and  NO}  in  this  energy 
range.  However,  Ranjan  and  Goodyear**  have  measured 
collisional  detachment  cross  sections  for  these  ions  in 
nitrogen  and  oxygen  at  energies  below  100  eV.  In  each 
case,  the  cross  section  values  are  well  below  our  higher 
energy  results,  but  are  still  increasing  at  100  eV.  At 
their  lower  energies,  there  appears  to  be  a  correlation 
between  the  magnitude  of  the  cross  section  and  the  neg¬ 
ative  ion  electron  attachment  energy. 

no;  is  known  to  exist  in  two  isomers  with  electron 
attachment  energies  of  about  2.4  and  3.9  eV.**  The  less 
stable  “excited” no;  is  formed  by  the  association  of 
an  O}'  with  a  NO  molecule  in  reactions  such  as  o;  +  NO 
->  0}  NO  *  O;,  and  the  more  stable  "ground  state”  NOf 
is  formed  In  reactions  such  as  NO^-f  NOf  NO;-t-NO. 

We  formed  NO;  ions  both  by  using  a  mixture  of  NO  and 
0}  and  by  using  pure  NO}  in  our  ion  source.  Mixtures 
of  all  three  gases  were  ^so  used.  Under  no  condition 
was  there  either  an  observable  departure  of  the  attenu¬ 
ation  curve  from  a  straight  line  or  a  dependence  of  the 
electron  loss  cross  section  on  the  gases  used  to  produce 
the  no;  .  Since  it  is  likely  that  the  NO;  produced  in  an 
NO/0}  mixture  would  be  primarily  the  less  stable  form 
and  that  the  NO;  pioduced  in  NO}  would  be  primarily 
the  ground  state  form,  the  electron  loss  cross  section 
appears  to  be  nearly  the  same  for  the  two  forms  of  NO; . 
This  is  consistent  with  the  apparent  lack  of  any  strong 
correlation  between  the  magnitude  of  the  loss  cross  sec¬ 
tion  and  the  negative  ion  electron  attachment  energy 


among  the  other  reactions  of  this  study.  However,  it  is 
also  possible  that  the  NO;  state  populations  were  nearly 
the  same  in  both  cases. 

Because  this  series  of  measurements  of  electron  loss 
cross  sections  was  peripheral  to  our  main  work  on  the 
merging  beams  apparatus  (ion-ion  mutual  neutraliza¬ 
tion),  the  scope  of  the  work  was  quite  limited.  Despite 
gaps  in  the  reactant  pairs,  as  Table  I  reveals,  the  re¬ 
sults  of  this  work  lend  weight  to  ideas  concerning  mech¬ 
anisms  of  collisional  detachment  of  electrons  from  neg¬ 
ative  ions.  In  all  reactions  studied  here,  the  cross  sec¬ 
tions  are  essentially  Independent  of  incident  ion  energies 
for  the  range  1-4  keV  (laboratory  energy)  except  for 
0‘+0}  for  which  charge  transfer  is  clearly  important. 
These  energies  are  well  above  threshold  for  detachment, 
but  the  incident  ion  speeds  are  always  less  than  2  x  10’ 
cm/sec.  This  is  below  the  range  of  validity  of  the  theo¬ 
ries**”*’  that  treat  the  collision  A'  +  B-  A  +  B  +  e‘  as  the 
scattering  of  a  loosely  bound  electron  moving  at  the 
speed  of  the  incident  ion.  At  the  Incident  ion  speeds  in 
our  studies,  we  would  expect  some  form  of  the  quasl- 
molecular  model**"**  to  apply.  In  this  model,  proposed 
lor  atomic  collisions  by  Bydin  and  Dukelskll,*'  the  po¬ 
tential  curves  for  the  ground  state  AB"  make  a  pseudo- 
crossing  with  the  ground  state  potential  curve  of  AB  at 
some  internuclear  separation  R,,  If  the  colliding  parti¬ 
cles  ^proach  closer  than  R,  the  electron  can  be  ejected 
with  some  probability  P  which  may  depend  on  the  inci¬ 
dent  ion  energy  and  the  potentials.  The  detachment 
cross  section  can  then  be  expressed  as 

Q  =  P(R.,W)ir/l*[l-K(Rj/H'],  (2) 

where  the  projectile  and  target  masses  are  A/,  and  Af}, 
respectively,  W  is  related  to  the  projectile  ion  kinetic 
energy  Tby  1^=  Af}r/(A/, +A/}),  and  K(R,)  is  the  poten¬ 
tial  energy  of  the  system  at  the  crossing  distance  R^. 

In  the  energy  range  of  our  experiments,  W  »  V'(R,) 
and  Eq,  (2)  reduces  to  Q=  PvR\ .  While  it  is  possible 
to  determine  P{W)  and  R,  for  some  systems, *•’**  the 


ENERGY  (keV) 


0.05  0.5  1  4  8  13 


BM= Bailey  and  Mahadevan.  *°  RG  =  Ranjan  and  Goodyear.  **  MS 
(charge  transfer  only)  =  Mathis  and  Snow.  ** 
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insensitivity  of  our  cross  sections  to  energy  suggests 
that  P  is  essentially  constant  in  our  energy  range.  Ap¬ 
parently  if  the  collision  energy  is  well  above  threshold, 
the  electron  is  ejected  with  unit  probability  when  the 
internuclear  separation  becomes  smaller  than  a  certain 
value.  The  cross  section  thus  shows  no  dependence  on 
energy.  This  model  will  hold  until  the  coUisional  veloc¬ 
ities  become  so  large  that  the  quaslmolecular  model  is 
no  longer  valid.  The  general  Increase  of  Q  with  pro¬ 
jectile  and  target  size  seen  in  Table  I  is  consistent  with 
this  model.  The  energy  dependence  of  Q  for  O'-t-Oi  is 
due  to  the  charge  exchange  channel.  This  channel  is 
also  possible  in  NO] ,  but  a  strong  energy  dependence 
of  Q  is  not  evident  there.  Finally,  there  appears  to  be 
no  systematic  dependence  on  negative  ion  electron  at¬ 
tachment  energy  for  the  cases  studied  here. 
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I .  Introduction 

The  two-body  ion-ion  mutual  neutralization  reaction 

A*+B“-A  +  B  (1) 

is  an  interesting  and  unusual  member  of  the  large  family  of  electron  capture  (charge 
transfer)  reactions.  It  is  of  theoretical  interest  because  it  has  certain  simplifying  char¬ 
acteristics  compared  to  most  ion-neutral  reactions.  It  is  of  practical  interest  be¬ 
cause  it  occurs  in  any  ionized  gas  where  negative  ions  are  formed  and  thus  has  im¬ 
portant  effects  in  electrical  discharges,  flames,  lasers,  and  in  the  earth’s  ionosphere. 

In  the  last  ten  years  considerable  effort  has  been  given  to  measurements  of  cross 
sections  and  rates  for  these  reactions,  and  to  improving  the  theoretical  methods  of 
treating  them.  This  paper  is  intended  to  review  the  results  of  the  various  experimen¬ 
tal  and  theoretical  methods  that  have  been  employed  to  date,  and  will  also  indicate 
some  problem  areas  that  are  not  yet  adequately  understood. 

The  basic  reaction  is  written  as  (1 ),  with  the  understanding  that  the  letters  A 
and  B  represent  either  atoms  or  molecules.  A  most  important  physical  characteristic 
of  the  reaction  is  that  the  two  reactant  ions  A^  and  B~  are  initially  under  the  in¬ 
fluence  of  the  Coulomb  potential,  whose  strong,  long  range  attractive  force  leads 
to  very  large  reaction  cross  sections  (exceeding  10~'^  cm^)  at  thermal  energies.  A 
second  important  characteristic  is  that  the  reactions  are  generally  exoergic  and  us¬ 
ually  leave  one  of  the  products  electronically  excited,  which  can  produce  radiation. 

If  the  excited  product  is  a  molecule,  it  may  be  in  an  unbound,  dissociating  state  and 
yield  other  fragment  products. 

Early  experimental  interest  in  two-body  ion-ion  neutralization  was  connected 
with  understanding  electron  and  ion  densities  in  electronegative  discharges  [  I  ]  and 
flames  [2] .  However,  the  major  efforts  both  theoretical  and  experimental  have  been 
stimulated  by  the  importance  of  the  reaction  in  the  D-region  of  the  earth’s  ionosphere 
(3-8) .  At  altitudes  of  60-90  km,  electrons  are  produced  in  the  daytime  by  photo¬ 
ionization.  These  electrons  may  be  removed  either  by  diss(x;iative  recombination 
such  as 
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e  +  N0"^-»N  +  0  (2) 

or  by  forming  a  negative  ion  in  attachment  reactions,  either  three  body 
e  +  0  +  M-0"+M  (3) 

or  dissociative 

e  +  O3  -» 0“  +  O2  .  (4) 

Following  attachment,  a  number  of  negative  ion-molecule  reactions  may  take 
place  to  change  the  nature  of  the  negative  ion  before  the  attached  electron  is  final¬ 
ly  returned  to  a  positive  ion.  During  the  daytime,  photodetachment  reactions  such 
as 

0~  +  /ii'-*0  +  e  (5) 

and  photodissociation  reactions  such  as 

CO3  +  hp  -*  CO2  +  0  (6) 

may  very  likely  occur  before  reactions (I)  removes  the  negative  ion,  starting  dif¬ 
ferent  parts  of  the  chain  of  reactions  over  again.  More  recently  it  has  been  found  [9| 
that  below  75  km  and  also  at  higher  altitudes  under  quiescent  conditions,  especially 
at  night,  both  the  positive  and  negative  ions  are  likely  to  undergo  hydration  reactions 
such  as 

NOj  +  HjO  +  M-NOj- H2O  +  M  .  (7) 

Subsequent  clustering  reactions  add  further  water  molecules,  and  other  reactions 
switch  the  original  parent  ions,  so  that  ions  eventually  evolve  into  combinations 
such  as  H30^  '  (H20)„  and  NO^  •  (H20)„,  whose  characteristics  are  far  removed 
from  those  of  the  original  ions.  On  the  other  hand,  during  disturbed  conditions, 
such  as  in  aurorae,  the  ions  tend  to  retain  a  simple,  unclustered  character  [10] . 

The  inverse  of  reaction  (1), 

A+B-A'^+B"  (8) 

cannot  easily  be  studied  because  the  neutral  products  of(l)  are  generally  left  in  ex¬ 
cited  states,  and  laboratory  measurements  of  (8)  must  usually  be  made  with  the  in¬ 
itial  neutral  species  A  and  B  in  their  ground  states.  The  production  of  H~  in  (ground 
state)  H  +  Mg  collisions  has  been  studied  at  collision  energies  above  5  keV  [1 1  ],  and 
cross  sections  for  O  +  Cs  -►  0  +  Cs^  have  been  measured  in  the  energy  range  1 80- 
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l600eV  [12].  Differential  cross  sections  for  Na+  I  -»Na^  +  1  have  been  measured 
at  energies  from  13  to  85  eV  (13) .  Several  studies  have  also  been  made  of  reactions 
such  as  A  +  BC  -►  AB*  +  C~,  but  these  are  even  farther  removed  from  the  true  in¬ 
verse  of(l). 

Theoretical  interest  in  the  ion-ion  neutralization  reaction  was  originally  due  to 
its  simplicity  among  heavy  particle  inelastic  collisions  that  can  be  understood  in 
terms  of  the  crossing  of  potential  energy  curves  1 14] .  The  curve-crossing  model  was 
first  proposed  in  1932  by  Landau  [  15] ,  Zener  Jib)  and  Stuckelberg  ]  17] ,  and  is 
commonly  referred  to  as  the  Landau -Zener  (LZ)  theory.  It  has  been  employed  ex¬ 
tensively  in  recent  years  in  theoretical  treatments  of  collisions  in  which  electronic 
transitions  occur.  The  simplifying  aspect  of  reaction  ( 1)  in  the  LZ  theory  stems  from 
the  fact  that  the  initial  state  of  the  reaction  (1),  pertaining  to  the  reactants  A^  +  B  , 
is  Coulombic  and  the  potential  energy  is  well  defined  at  large  intcrnucicar  separa¬ 
tions/?,  by  L(/?)=  -e^lR. 

In  zeroth-order  approximation,  ignoring  interactions  between  electronic  states, 
this  Coulomb  potential  is  intersected  by  a  number  of  potentials  of  excited  covalent 
states  of  the  molecule  AB,  which  arc  essentially  flat.  Thus  the  positions  of  the  cross¬ 
ings  and  the  slopes  of  the  zeroth-order  curves  at  the  intersections  are  easily  ob¬ 
tained.  These  quantities  enter  directly  into  the  Landau-Zener  ( LZ)  theory,  and 
the  simplicity  in  obtaining  them  for  reactions  of  type  ( 1 )  is  attractive  from  the 
theorists’  standpoint.  The  various  theoretical  treatments  will  be  described  in  detail 
later  in  sections  3  and  5. 

Ion-ion  mutual  neutralization  has  been  reviewed  briefly  by  Sayers  [18]  and 
Bates  1 19] .  In  a  recent  article  in  this  series  [20] ,  Bates  discusses  the  general  topic 
of  recombination,  with  emphasis  on  the  historical  development,  applications  and 
theory.  Mahan  [21  ]  has  also  treated  the  general  case  of  recombination  of  gaseous 
ions,  including  three-body  (“Thomson”)  recombination  A*  +  B~  +  M  -*  (A  +  B)  +  M. 
Flannery  [22]  has  made  a  definitive  review  of  the  theory  of  three-body  ion-ion 
recombination.  For  a  treatment  of  all  of  these  and  related  areas  the  reader  is  refer¬ 
red  to  the  comprehensive  books  Collision  Phenomena  in  Ionized  Gases  by  McDaniel 

[23]  and  Electronic  and  Ionic  Impact  Phenomena  by  Massey,  Burhop  and  Gilbody 

[24] . 

In  this  article  we  deal  only  with  the  two-body  reactions,  with  the  goal  of  present¬ 
ing  a  reasonably  thorough  picture  of  it  as  a  two-body  inelastic  charge  transfer  reac¬ 
tion.  We  shall  review  the  various  expcrimen.al  and  theoretical  studies  that  have  been 
made,  in  an  attempt  to  understand  the  physics  of  the  reactions  as  well  as  to  obtain 
reliable  rate  constants  that  are  required  for  various  applications  in  ionized  gases. 


2.  Experimental  techniques 

2.1.  Measurement  of  ion  density  decay 

The  first  approaches  to  ion-ion  recombination  studies  involved  determinations 
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of  the  rates  of  decay  of  ion  densities  in  reaction  chambers  following  some  form  of 
ionization  process.  In  an  ideal  case,  the  ion  densities  are  monitored  for  a  period  of 
time  after  all  of  the  initially  free  electrons  in  the  volume  have  been  removed  by  at¬ 
tachment  or  recombination,  and  under  conditions  such  that  the  spatial  distributions 
of  the  ions  are  uniform  and  no  other  loss  mechanism  exists  except  ion-ion  neutrali¬ 
zation  (either  two  or  three  body).  Also,  ideally  only  one  molecular  or  atomic  form 
of  either  positive  or  negative  ion  is  present  in  the  volume,  so  that  the  ion  identities 
are  known  and  only  one  reaction  is  involved.  Under  these  conditions  space  charge 
neutrality  requires  that  the  positive  and  negative  ion  densities  p'*^  and  p  are  equal. 
Then  the  rate  of  change  of  either  is  given  by 

dp/ (it  =  -  ap^  ,  (9) 

and  the  ion  densities  themselves  have  time  dependences  expressed  by  p“  ’  = 

Pjj '  +  at,  where  Pq  is  the  density  at  r  =  0  when  the  measurement  is  started.  At  any 
pressure,  the  reaction  rate  coefficient  a  can  be  obtained  from  the  slope  of  a  plot  of 
p  *  versus  t.  In  the  pressure  range  where  both  two-body  and  three-body  processes 
are  important,  the  effective  neutralization  rate  or  may  be  expressed  as  a  =  Oj  +  03, 
the  sum  of  the  respective  two-body  and  thre6-body  rate  coefficients.  Three-body  re¬ 
combination  can  be  viewed  as  the  formation  of  an  unstable  complex  (A*  +B  -*AB*) 
which  is  then  stabilized  into  AB  or  A  B  by  collision  with  a  gas  molecule.  Thus  the 
three-body  rate  coefficient  03  will  be  linear  in  pressure  until  the  pressure  becomes 
high  enough  that  the  time  between  collisions  approaches  the  lifetime  of  the  com¬ 
plex  AB*,  when  03  will  lose  its  linearity  with  pressure.  The  two-body  rale  03  can  in 
principle  be  determined  from  the  zero  pressure  intercept  when  a  is  plotted  versus 
pressure,  but  this  will  be  valid  only  when  it  is  established  that  03  is  linear  with  pres¬ 
sure.  In  many  cases,  even  when  03  is  linear  with  pressure,  Oi  may  be  much  smaller 
than  03,  and  the  extrapolation  method  of  determining  03  is  not  very  accurate. 

Ion  densities  have  been  monitored  by  several  techniques.  Some  of  the  earliest 
studies,  by  Yeung  [  l|  and  Greaves  |25| ,  utilized  a  measurement  of  the  dielectric 
constant  of  the  afterglow  plasma  following  ionization  by  r.f  fields.  Shifts  in  the 
resonant  frequencies  of  r.f  sampling  probes  were  the  measured  parameters.  In  a  re¬ 
lated  method,  Knewstubb  and  Sugden  (2)  measured  the  change  in  the  Q  of  a  re¬ 
sonant  r.f  circuit,  caused  by  ionization  in  an  alkali  water-vapor  flame.  They  as¬ 
sociated  this  shift  with  the  electrical  conductivity  and  thus  determined  the  ion 
densities.  Recently,  a  more  sophisticated  experiment  by  Eisner  and  Hirsh  [26|  de¬ 
termined  the  conductivity  of  air-like  mixtures  following  ionization  by  high  energy 
electrons. 

Most  of  these  ion  density  decay  measurements  have  been  made  using  apparatures 
in  which  ion  currents  from  the  afterglow  plasmas  have  been  collected,  usually  using 
weak  electric  fields  to  draw  the  ions  to  collector  plates,  and  have  only  determined 
three-body  "volume”  recombination  coefficients  at  pressures  above  10  torr.  How¬ 
ever,  Mahan  and  coworkers  [27-29]  made  extensive  measurements  at  pressures  down 
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to  less  than  1 0  torr,  which  allowed  them  to  deduce  two  body  recombination  coef¬ 
ficients. 

In  practice,  the  ideal  conditions  listed  at  the  start  of  this  section  never  exist  and 
the  analysis  required  for  proper  interpretation  of  the  results  is  very  complicated. 
Consequently,  much  of  the  early  work  was  improperly  analyzed.  In  his  recent  review, 
Mahan  (21  ]  has  summarized  the  complexities  of  dealing  with  real  systems.  Only  a 
brief  discussion  of  some  of  these  complications  will  be  given  here. 

Diffusion  losses  modify  the  loss  rate  equation  (9)  so  that  it  becomes 

3p/ar  =  -  0^2  _  /)v2p  (10) 

where  D  is  the  effective  diffusion  coefficient  which  is  often  assumed  equal  for  po¬ 
sitive  and  negative  ions.  McGowan  [30]  analyzed  diffusion  effects  for  a  plane  paral¬ 
lel  ionization  chamber  and  applied  the  results  to  his  measurements  [31 1.  Fisk,  Ma¬ 
han  and  Parks  [29]  used  a  numerical  solution  to  account  for  their  diffusion  effects. 

In  addition  to  these  diffusion  losses,  which  can  be  accounted  for,  inhomogeneities 
in  the  initial  ion  densities  can  lead  to  errors  in  the  interpretation  of  the  decay  rates. 

A  common  problem  in  most  experiments  is  that  more  than  one  ion  species  of 
each  charge  exists  in  the  afterglow.  In  many  gas  mixtures,  several  ion-molecule 
and  charge  transfer  reactions  can  take  place  before  the  most  stable  species  of  each 
is  reached.  Until  this  stabilization  is  reached,  ion-ion  neutralization  occurs  between 
several  species,  and  the  effective  rate  is  some  average  of  the  individual  rates.  Since 
the  diffusion  coefficients  may  be  different  for  each  species,  the  analysis  becomes 
difficult.  Mass  analysis  has  been  used  to  sample  the  ions  in  the  more  reliable  ex¬ 
periments.  Thus,  Greaves  [25|  was  able  to  identify  I2  +  I~  as  the  only  reacting  pair 
in  his  chi'mber.  Similarly,  Fisk,  Mahan  and  Parks  [29]  identified  their  reactants  as 
+1  ,  Tl*  +  Cr,  Tl*  +  Br  only,  in  pure  vapors  of  Til,  TICI  and  TIBr,  respectively, 
at  low  pressure,  while  at  higher  pressures  (and  temperatures),  the  dominant  ions  be¬ 
come  TI2I*  and  Tllj,  etc.  On  the  other  hand,  mixtures  of  Til  and  NO2  produced  a 
complicated  spectrum  whose  effective  recombination  coefficients  were  not  easily 
interpretable.  Mass  sampling  provides  a  clear  interpretation  when  only  one  reaction 
is  controlling  the  ion  decay,  but  the  results  become  less  useful  when  a  variety  of 
ions  is  found. 

In  examining  air-like  mixtures  that  contained  several  ions  of  each  charge,  whose 
relative  densities  depend  on  the  length  of  time  since  irradiation  by  I  MeV  electrons, 
Hirsh  and  Eisner  [26,32|  used  a  theoretical  model  to  obtain  rate  coefficients  from 
their  mass  spectrometer  data.  Ions  were  sampled  from  air-like  mixtures  at  2-22  torr 
following  bombardment  by  I  MeV  electrons.  A  varying  number  of  ions  of  several 
species  were  found,  whose  relative  popi’lations  depended  on  the  gas  mixture,  length 
of  irradiation  time,  and  pressure.  The  theoretical  model  was  used  to  deduce  two- 
body  neutralization  rates  for  NO*  +  NO2  and  NO*  +  NOJ  from  the  sampled  ion 
current  data.  Total  ion  density  decay  rates  were  determined  from  plasma  conduc¬ 
tivity  measurements  using  an  r.f  probe.  This  technique  is  particularly  appealing  for 
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studies  of  atmospheric  ion  species  because  it  can  approximate  some  of  the  cond; 
tions  that  must  exist  for  the  natural  formation  of  these  ions.  However,  proper  data 
analysis  requires  a  model  that  must  account  for  different  diffusion  rates,  wall  ef¬ 
fects,  and  possibly  competing  reactions.  An  accurate  model  can  thus  become  veiy 
complicated,  and  tests  of  its  validity  are  difficult.  Since  the  early  measurements  of 
Eisner  and  Hirsh  (26| ,  the  model  has  been  further  advanced  [33| . 

2.2.  Merged  beams 

In  order  to  overcome  some  of  the  uncertainties  inherent  in  most  bulk  measure¬ 
ments  of  ion  densities  described  in  the  previous  sections,  the  mc^i  corr^ticated 
“merged  beams"  method  was  developed  for  ion-ion  reactions  at  Stan*'  rj  Research 
Institute  [34, 3S| .  In  this  method,  two  mass-selected  ion  beams  moving  in  the  same 
direction  are  superimposed  over  a  known  path  length,  and  then  separated.  The  ions 
each  travel  at  high  speeds  in  .lie  laboratory  frame  of  reference,  with  kinetic  energies 
of  several  keV,  but  the  relative  speeds  in  the  center  of  mass  frame  can  be  made  very 
small,  j’pproaching  thermal  velocities.  The  identities  of  the  ions  can  be  established, 
and  the  relative  collision  energies  can  be  accurately  controlled  over  a  wide  range. 
Thus,  the  cross  section  energy  dependence  as  wi?ll  as  its  absolute  magnitude  can  be 
established,  providing  more  insight  into  the  details  of  the  electron  transfer  reactions 
than  is  possible  by  bulk  measurements. 

The  merged  beam  technique  was  first  used  to  study  ion-neutral  reactions  and 
has  been  reviewed  by  Neynaber  (36) .  In  addition  to  the  advantage  of  mass  selection, 
the  technique  offers  a  very  precise  determination  of  the  relative  kinetic  energies  of 
the  colliding  ions,  as  may  be  seen  from  the  following  brief  analysis. 

The  relative  energy  Ej  of  ions  in  two  beams  traveling  at  speeds  Uj  and  1*2  in  the 
laboratory  and  at  an  an^e  6  with  respect  to  each  other  is 

£’,  =  i  i  j  +  uj  2u,  1^2  cos  0) .  (11) 

where  n  is  the  reduced  mass.  In  a  merged  beams  apparatus,  ideally  d  =  0  and  the 
relative  energy  £',o  is  given  by 

M  (1^1  -  *^2)^  =  )’  -  (^2/^2)’  02) 

where  £,•  and  Af,  refer  to  the  laboratory  energy  and  mass  of  each  ion.  The  precision 
in  the  relative  energies  results  from  a  "deamplification”  in  transforming  the  laboratory 
energies  into  the  relative  (center  of  mass)  collision  energies.  This  effect  is  easily  seen 
in  the  case  of  equal  masses  rn|  =  m2  and  small  energy  differences  Af  =  £2  ^  ^1-  )n 
this  case  eq.  (12)  reduces  to  £,  ^  (A£)2/8f,  where  £  is  the  mean  lab.  energy.  A  simi¬ 
lar  deamplification  in  the  general  case  reduces  the  normal  energy  spreads  of  a  few  eV 
in  each  of  the  beams  to  very  small  amounts  in  the  center-of-mass  system.  In  fact,  in 
a  practical  experiment,  when  the  average  relative  energy  is  in  'he  range  below  1  eV. 
the  greatest  uncertainty  in  the  relative  energies  is  due  to  imperfect  alignment  and  col- 
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limation  of  the  beams  [35, 37] .  This  effect  can  be  understood  with  the  aid  of  eq. 
(11),  where  0  represents  the  angle  between  any  two  colliding  ions.  When  0  is  small, 
eq.  (11)  becomes 

fjQ  is  the  relative  energy  of  the  two  beams  determined  from  the  difference  in  the 
laboratory  energies,  as  one  would  obtain  from  eq.  ( 1 2).  Thus,  because  of  the  trans¬ 
verse  velocity  components,  the  average  energy  is  increased,  and  an  uncertainty  in  0 
produces  an  uncertainty  in  E^.  Practical  lower  limits  to  E^  are  generally  near  0.1  eV. 

The  original  apparatus  at  Stanford  Research  Institute  (SRI)  [34,35]  has  under¬ 
gone  several  mternal  modifications,  and  is  presently  in  a  form  represented  schema¬ 
tically  in  fig.  1 .  The  two  beams  are  generated  in  duoplasmatron  ion  sources,  mass- 
selected  in  Wien-type  EX  B velocity  filters,  focused,  and  merged  in  a  common 
magnet.  The  superimposed  beams  then  enter  an  ultra-high  vacuum  chamber  and 
flow  together  over  the  interaction  path  until  they  are  separated  by  electrostatic 
deflection  (at  “Demerger  B”  in  fig.  1)  and  finally  collected  for  current  measure¬ 
ment.  The  fast  neutral  particles  formed  from  the  beams  along  the  interaction  path 
continue  undeflected  until  they  are  stopped  and  detected  by  secondary  electron 
emission.  Complete  collection  is  possible  because  all  reaction  products  are  in  a  nar¬ 
row  cone  surrounding  the  beam  axis. 

The  pressure  m  the  interaction  chamber  is  maintained  below  2  X  10'^  torr 
during  the  experiments  by  titanium  sublimation  pumping,  backed  by  an  oil  diffu¬ 
sion  pump  (the  ultimate  pressure  with  beams  off  is  about  5  X  10'**^  torr).  How¬ 
ever,  even  at  2  X  10"^  torr,  the  fast  neutrals  in  the  beam  produced  by  charge  trans- 


Fig.  I .  Schematic  of  the  merged  beams  apparatus  at  Stanford  Research  Institute. 
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fer  (A*  +  G  -+  A  +  G"*^)  and  collisional  detachment  (B~  +  G-*^B  +  G  +  e)  reactions 
with  the  background  gas  (G)  outnumber  the  ion— ion  neutralization  products  by  one 
or  two  orders  of  magnitude.  In  order  to  eliminate  these  beam -gas  contributions 
from  the  neutral  product  signal,  the  two  beams  are  chopped  at  different  rates  and 
the  signal  is  coherently  detected  at  the  difference  frequency  [3S| . 

To  define  the  interaction  path  length,  and  to  eliminate  the  interference  signal 
which  arises  from  beam-beam  interactions  as  the  beams  are  being  merged  in  the 
magnet  chamber,  the  beams  may  also  be  separated  by  a  pair  of  platej  (“Demerger 
A”)  located  just  inside  the  interaction  chamber.  The  effective  interaction  path  length 
is  thus  the  distance  between  the  demergers  (33  cm). 

Changes  in  relative  energies  are  made  in  two  ways.  Small  changes  over  a  limited 
range  are  easily  made  by  changing  the  potential  on  an  electrostatic  cage  which  sur¬ 
rounds  the  merged  beams  over  the  path  between  the  two  sets  of  deflection  plates. 

A  pair  of  high  transparency  plane-parallel  grids,  oriented  perpendicularly  to  the 
beam  path,  is  located  Just  in  front  of  the  cage.  One  grid  is  electrically  attached  to 
t'le  electrostatic  cage  and  forms  the  entrace  “window”  to  the  cage.  The  other  is 
ground  and  is  located  about  0.5  cm  ahead  of  the  cage.  Thus,  only  the  uniform 
field  between  the  grids  is  experienced  by  the  ions  as  they  enter  the  cage.  A  poten¬ 
tial  on  the  cag'.'  will  change  the  relative  laboratory  energies  by  21^^,.  In  practice, 
potentials  of  up  to  ±  400  V  can  generally  be  used  without  seriously  affecting  the 
beam  trajectories,  allowing  a  relative  energy  range  of  order  of  10-50eV  to  be  easily 
covered.  Other  sets  of  energies  can  be  covered  by  changing  one  of  the  main  beam 
energies,  which  requires  changing  its  entrance  angle  at  the  magnet  in  order  that  it 
have  the  same  exit  angle  as  the  other  beam.  This  may  be  done  either  (35)  by  chang¬ 
ing  the  angle  of  the  appropriate  ion  source  chamber  (both  are  rotatable  about  the 
center  of  the  magnet),  or  by  using  two  sets  of  deflectors  [38]  located  in  front  of 
the  merging  magnet,  labeled  “Deflectors”  I  and  2  in  fig.  1 .  Once  the  beam’s  energy 
is  cha.nged  and  the  beam  is  realigned,  the  electrostatic  cage  potential  may  be  varied 
to  cover  a  new,  partially  overlapping,  range  of  relative  energies.  The  ability  to  cover 
the  sam?  relative  energy  with  more  than  one  laboratory  beam  energy  and  entrance 
angle  affords  a  check  on  the  reliability  of  the  data. 

Proper  beam  focusing  and  adjustment  require  great  care.  The  neutral  product 
signal  depends  (35)  on  the  beam  cunent  density  overlap  integral  where  the 

J's  represent  the  current  densities  and  the  integral  is  taken  over  all  the  volume  V 
common  to  the  two  beams  along  the  interaction  path.  Focusing,  adjustment,  and 
an  approximate  solution  to  the  overlap  integral  are  facilitated  with  the  use  of  two 
movable  beam  flags  containing  small  apertures,  each  0.3  Ciis  in  diameter.  These  can 
be  moved  accurately  into  position  centered  on  the  beam  axis  at  each  end  of  the 
interaction  path  next  to  the  larger,  fixed  1.3  cm  diameter  apertures  at  each  end,  as 
described  in  [3S|  and  [39]. 

The  signal  also  depends  on  the  average  secondary  electron  emission  coefficients 
7  appropriate  to  the  fast  neutral  products  that  strike  the  neutral  product  detector. 
This  detector  is  a  stainless  steel  plate,  oriented  at  45°  to  the  beam  direction  in  order 
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to  increase  the  secondary  electron  yield.  Hie  secondary  electron  currents  that  flow 
from  this  plate  to  the  walls  of  the  surrounding  detector  enclosure,  which  are  biased 
at  a  positive  voltage,  provide  the  neutral  product  signal.  The  value  assigned  to  j  is 
the  average  of  y*  and  y~,  the  coefficients  corresponding  to  the  two  ion  beams, 
which  are  measured  frequently  during  the  experiments.  This  method  of  determining 
7,  and  its  reliability,  are  discussed  in  !39] .  It  is  regarded  to  be  accurate  to  within 
10%  of  the  actual  effective  y  of  the  neutral  products,  direct  measurement  of  which 
is  experimentally  inaccessible. 

In  order  to  reach  minimum  energies,  the  average  laboratory  speeds  of  the  two 
beams  must  be  made  equal.  Under  this  condition  the  momenta  of  the  positive  and 
negative  ions  are  proportional  to  the  respective  masses,  and  the  angles  through 
which  the  beams  are  deflected  in  the  common  merging  magnet  are  roughly  inversely 
proportional  to  the  ion  masses.  Thus,  the  heavier  ion  is  deflected  less  than  the 
lighter  ion,  and  correspondingly  less  mass  (momentum)  selection  is  afforded  the 
heavier  ion.  In  most  cases  several  ions  of  the  same  charge  may  exist  in  the  beam 
produced  by  the  ion  source,  and  often  the  merging  magnet  does  not  provide  ade¬ 
quate  mass  selection.  As  a  result,  Wein  £X  B  velocity  filters  were  installed  along 
each  beam  path  to  provide  primary  mass  selection. 

Although  the  average  laboratory  speeds  in  the  two  beams  may  be  equal,  the 
lower  limit  to  the  effective  relative  energy  in  the  beams  is  non-zero,  due  to  velocity 
components  transverse  to  the  beam  direction  stemming  from  incomplete  collima- 
tion  and  imperfect  focusing  (see  eq.  ( 13)).  In  the  SRI  apparatus,  the  lower  limit  to 
the  center-of-mass  energy  is  in  the  range  0.1  -0.2  eV  due  to  these  effects. 

The  total  energy  range  covered  in  most  measurements  extends  from  about  0. 1 5 
eV  to  about  200  eV.  This  large  dynamic  range  has  aided  the  development  and  test¬ 
ing  of  fairly  effective  but  still  uncomplicated  theoretical  models  for  both  atomic 
ions  [40,4 1 1  and  molecules  [42 j .  However,  more  accurate  theories  will  require 
further  refined  measurements  of  final  states  since  only  total  cross  sections  have 
been  subject  to  extensive  study. 

In  order  to  predict  thermal  reaction  rates,  which  are  required  for  applications 
in  aeronomy,  the  SRI  group  has  used  an  extrapolation  technique  to  extend  their 
measured  cross  sections  to  lower  energies  [37-39] .  This  technique  will  be  discuss¬ 
ed  in  section  3. 

A  different  configuration  of  a  merged  beams  apparatus  has  been  used  by  Weiner, 
Peatman  and  Berry  [4.^, 44]  at  the  University  of  Chicago  to  study  optical  emissions 
from  the  final  states  jf  Na*  +  0~  neutralization  reactions.  Figure  2  is  a  schematic 
drawing  of  that  apparatus.  Here  each  ion  beam  is  deflected  through  180°  in  the 
merging  magnet.  Photons  can  be  observed  th.'<'ugh  two  windows,  one  on  the  top 
and  one  on  the  side  of  the  vacuum  chamber.  Interference  filters  are  used  to  select 
the  desired  emission  line,  and  the  light  is  focused  by  a  lens  onto  a  photomultiplier. 
Individual  photons  are  counted  and  a  sequence  of  beam-chopping  and  count-gating 
pulses  permits  various  background  counts  to  be  subtracted  from  the  total  photon 
counts.  The  signal  to  noise  ratio  is  about  0.1. 
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Fig.  2.  Schematic  of  the  merged  beams  apparatus  of  Weiner,  Peatman  and  Berry  |43,44|. 

The  primary  difnculty  in  obtaining  absolute  neutralization  cross  sections  from 
this  experiment  is  calibration  of  the  sensitivity  of  the  photon  detection  system.  The 
optical  collection  efficiency  is  estimated  geometrically,  and  the  single  photon  count¬ 
ing  efftciency  of  the  photomultiplier  is  assumed  to  be  the  same  as  the  quantum  ef¬ 
ficiency  quoted  by  the  manufacturer.  It  has  recently  been  shown  [45]  that  these 
quantities  can  differ  substantially.  Therefore  the  uncertainty  in  the  photon  count¬ 
ing  sensitivity  is  difficult  to  assess,  and  the  absolute  value  of  the  cross  sections  de¬ 
termined  by  this  experiment  are  subject  to  some  question. 

Since  this  experiment  observes  optical  emission  from  the  neutral  products  of 
an  ion-ion  neutralization  reaction,  it  provides  direct  information  about  the  final 
states  of  these  reactions  and  the  energy  dependences  of  the  various  reaction  chan¬ 
nels.  Thus  important  new  information  can  be  made  available  for  the  theoretical 
understanding  of  the  collision  processes,  and,  as  will  be  discussed  in  section  4,  some 
surprising  results  have  been  obtained. 

2.3.  Inclined  beams 

The  first  measurements  of  one  of  the  most  theoretically  interesting  ion-ion 
neutralization  reactions,  -♦  H  +  H,  were  performed  by  Rundel,  Aitken 

and  Harrison  [46]  and  by  Gailey  and  Harrison  [47]  at  the  United  Kingdom  Atom 
tc  Energy  Authority  Laboratory  in  Culham,  England,  using  an  “inclined  beams” 
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Fig.  3.  Schematic  of  the  inclined  beams  apparatus  of  Harrison  and  co-workers  |46  48|. 


technique.  A  schematic  of  their  apparatus  is  shown  in  fig.  3.  The  positive  and  ne¬ 
gative  beams  intersect  each  other  at  an  included  angle  0  (see  eq.  ( 1 1 ))  of  20°.  Fol¬ 
lowing  the  collision  region,  the  proton  beam  is  further  deflected  into  a  Faraday  cup. 
while  neutrals  formed  proton  interactions  with  the  H”  beam  and  the  background 
gas  continue  into  the  neutral  detector.  Chopping  of  both  beams  is  employed  to  al¬ 
low  subtraction  of  background  effects. 

While  this  technique  does  not  permit  measurements  at  nearly  thermal  relative 
energies,  as  in  the  merged-beams  methods  when  0  =  0°,  it  does  allow  the  energy 
range  to  be  extended  to  much  higher  energies.  Thus  the  +  H~  studies  were  car¬ 
ried  out  at  large  center  of  mass  energies,  between  1 25  and  5000  eV.  At  these  re¬ 
latively  high  energies  it  is  important  to  separate  the  products  from  the 
collision  since  electron  detachment,  H*  +  H~  H*  +  H  +  e,  becomes  increasingly 
likely.  By  observing  only  neutrals  formed  from  the  positive  beam,  this  experiment 
measures  li  e  ion-ion  mutual  neutralization  cross  section  without  possible  interfer¬ 
ence  from  electron  detachment.  This  apparatus  was  also  used  to  make  measure¬ 
ments  on  He*  +  H~  [48| ,  which  was  also  studied  at  SRI  [40] . 

2.4.  Discussion  of  experimental  lechniques 

The  three  basic  techniques  that  have  been  used  to  study  ion-ion  neutrali7.ation 


40 


14  J.T.  Moseley,  R.H.  Olson  and  J.R.  Peterson.  Ion  ion  mutual  neutralization 

are  complementary.  The  bulk  gas  techniques  allow  measurements  to  be  made  at 
thermal  energy  between  ions  that  are  in  most  cases  likely  to  be  in  their  ground 
states.  In  order  to  unambiguously  determine  the  reaction  rate  for  a  particular  pair 
of  ions,  however,  it  is  generally  necessary  to  use  a  mass  analyzer.  Even  then,  analysis 
becomes  very  complex  if  more  than  two  species  of  ions  are  present. 

The  merged  beams  technique  allows  measurements  on  well  identified  pairs  of 
ions  over  a  wide  energy  range  (0.1  to  several  hundred  eV).  Absolute  cross  sections 
can  be  measured  as  a  function  of  energy,  and  by  observing  optical  emission  the 
final  neutral  states  can  be  investigated.  Absolute  accuracy  is  so  fa'  limited  to  about 
30%,  mainly  by  the  difficulty  in  determining  the  overlaps  of  the  ion  beams,  and 
is  even  further  restricted  when  photon  measurements  are  used.  Further,  for  mole¬ 
cular  ions,  uncertainty  in  the  vibration."!  and  electronic  states  of  the  interacting 
ions  leads  to  some  ambiguity  in  the  results,  as  will  be  discussed  in  more  detail  later. 

The  inclined  beams  techniques  (as  well  as  any  other  “crossed  beam”  measure¬ 
ment)  has  similar  advantages  and  disadvantages  to  the  merged  beams  technique.  It 
cannot  approach  the  low  energy  limits  of  the  merged  beams  but  can  distinguish  be¬ 
tween  neutral  products  from  the  positive  and  negative  beams. 

3.  Curve  crossing  model  for  ion-ion  mutual  neutralization 
3.1.  Introduction 

A  theoretical  model  that  has  been  used  with  some  success  to  calculate  the  ion- 
ion  mutual  neutralization  total  cross  section  is  the  Landau-Zener  [15-17]  curve 
crossing  theory.  A  general  description  of  the  application  of  this  model  to  the  ion- 
ion  neutralization  problem  was  first  given  by  Bates  and  Massey  [14]  and  subse¬ 
quently  reviewed  by  Massey  [49] .  Later,  Magee  [50]  applied  the  Landau-Zener 
theory  to  the  multicrossing  case  ofO*  +  M”  ->  0  +  M  where  M  was  O  and  02- 

Tliis  theoretical  approach  can  be  discussed  with  the  aid  of  fig.  4.  For  this  dis¬ 
cussion  it  will  be  assumed  that  both  ions  are  atomic.  The  initial  state  of  the  system 
is(A^  +  B~)  at  internuclear  separation/?  =<».  If  the  ions  approached  each  other 
without  interacting  except  through  the  Coulomb  force,  the  potential  energy  of  the 
system  would  follow  the  Coulomb  potential  curve  -c^jR.  In  fig.  4  this  Coulomb 
curve  is  “crossed”  by  those  of  several  excited  neutral  states.  These  neutral  state 
curves  approach  /?  =  0  horizontally  until  A  and  B  begin  to  interact  at  close  separa¬ 
tions.  usually  about  4  to  6  atomic  units.  Consider  the  crossing  labe'ed  “2”  in  fig. 

4.  As  (A*  +  B~)  approaches  such  a  crossing  with  a  neutral  state  (.\  +  B)  of  the  same 
symmetry,  an  interaction  between  the  states  can  occur,  which  causes  an  actual  cros¬ 
sing  to  be  avoided,  as  indicated.  During  the  collision,  as  R  decreases  past  the  cros¬ 
sing  distance  /?^,  the  system  has  some  probability  p  of  making  a  transition  to  the 
other  state.  This  same  probability  exists  as  the  collision  partners  separate.  The  theo¬ 
ries  of  Landau  [  15] ,  Zener  1 16]  and  Stueckelberg  [  17]  provide  a  way  of  calculat¬ 
ing  this  probability,  which  will  be  discussed  later  in  this  section.  For  a  single  cros¬ 
sing  such  as  “2,”  the  probability  Z’  of  a  transition  from  (A^  +B~)  to  (A  +  B)  is/’  = 
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rif!.  4.  Simplified  potential  energy  eurves  for  a  “typical”  atomic  ion  system. 


2p(l-p).  and  the  result  for  multiple  crossings  follows  logically. 

The  ion-ion  neutralization  reaction  is  appealing  theoretically  because  it  offers 
the  simplest  test  of  the  LZ  theory.  The  crossing  points  and  the  difference  in  slopes 
of  the  two  non-interacting  curves  at  the  crossing,  which  arc  needed  to  calculate  the 
transition  probability,  can  be  accurately  determined  since  the  Coulomb  curve  is 
known  and  the  neutral  potentials  are  essentially  flat. 

3.2.  Landau -Zener  theory 

The  transition  probability  for  a  crossing  between  an  initial  state  i  and  a  final  state 
f  is  given  by  the  Landau-Zener  model  (15  17)  as* 


Pif  =  exp(-2u^/uft)  (14) 

where  Vf,  is  the  radial  velocity  for  impact  parameter  h  at  the  crossing  point  R^,  and 

u,  =  jr//jf2(/?J/|C;-f^fl.  (15) 

In  eq.  (15)  and  Vf  are  the  initial  and  final  state  potentials,  respectively,  //jf  is 
the  coupling  matrix  element  between  the  states,  and  is  dVJdR  evaluated  at  R^. 
The  total  cross  section  is  given  by 

Qif=2nfl^bPif{b)db.  (16) 

where  the  impact  parameter  b  is  related  to  the  angular  momentum  /  and  the  wave 
number  .'c  fy 

*  Atomic  units  will  be  used  throughout  this  section  and  section  S. 
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/>  =  (!  + i)/*.  (17) 

The  impact  parameter  for  the  crossing  point  is  given  by 
I 

b^=R^[\-V,(R^)IEy  (18) 

In  the  special  case  of  only  one  crossing,  Pd  has  the  form,  as  mentioned, 

Pir-W-Pif)-  09) 

In  a  multicrossing  system,  however,  allowance  must  be  made  for  the  change  in  P^f 
due  to  other  curve  crossings  [SO]. 

For  the  case  at  hand,  reaction  ( I ), 

+  A  +  B  +  AA’,  (1) 

where  A  or  B  may  each  or  both  be  in  excited  electronic  states,  the  initial  potential 
curve  is  Coulomb,  and  the  final  state  is  that  of  separated  neutral  atoms.  For  mod¬ 
erate  energies,  important  crossings  occur  only  over  a  limited  range  of  internuclear 
separation  /?.  At  large  distances,  the  coupling  matrix  elements  die  off  exponential¬ 
ly,  while  for  close  encounters,  transitions  occur  only  at  high  velocities  and  small 
impact  parameters,  and  these  make  only  a  small  contribution  to  the  total  cross  sec¬ 
tion.  Favored  aossings  occur  from  about  10  to  50  a.u.  For  these  separations,  to  a 
good  approximation  the  interactions  between  the  neutral  atoms  may  be  neglected 
compared  to  the  Coulombic  one.  Then 

AE=-V^(R^)  =  R;'  (20) 

and 

\Vy-V'(\=R-'^ .  (21) 

The  value  of  A£'  can  be  calculated  directly  from  the  ionization  potential  of  A, 
I.P.(A).  the  electron  affinity  of  B,  E.A.(B),  and  the  levels  of  excitation  of  the  neu¬ 
tral  products,  E{\),  and  ffB).  For  ground  state  products  M'  is  given  by 

A£=  [l.P.(A)-E.A.(B)|  .  (22) 

In  more  defined  treatments  [51-53]  of  reaction  (1),  the  polarizabilities  of  the 
initial  and  final  states  are  taken  into  account  to  obtain  more  accurate  evaluations 
of  eqs.  (20)  and  (21 ).  This  effect  becomes  very  important  at  small  crossing  distances, 
but  at  the  separations  of  importance  here  (10-50  a.u.)  eqs.  (20)  and  (2 1 )  are  adequate. 
For  the  alkali-oxygen  systems,  this  has  been  verified  by  van  den  Bos  [53]  and  for 
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the  H'''  -t-  H~  reaction,  the  values  of  Bates  and  lutwis  [51]  do  not  deviate  signifi¬ 
cantly  from  the  predictions  of  the  above  equations. 

Therefore,  from  eqs.  (IS)  and  (21),  we  obtain 

Fx=7r/?2//.f2(/?x)  (23) 

and  from  (18) 

(24) 

The  radial  velocity  at  is  given  by 

Ui=u[l+(f?xf)-'-(6//?x)^l^  .  (25) 

where  v  is  the  incident  velocity.  From  a  knowledge  of  the  electron  affinity  of  A,  the 
ionization  potential  of  B,  and  the  electronic  energy  levels  of  A  and  B,  the  AE's  of 
reaction  (1)  for  transitions  to  the  various  states  may  be  calculated  by  eq.  (22).  Then, 
by  eq.  (20),  the  crossing  point  for  each  AE  is  known.  Now  all  the  quantities  neces¬ 
sary  for  the  calculation  are  realized,  except  for  the  coupling  matrix  elements. 

These  matrix  elements  are  difficult  quantities  to  obtain  by  ab  initio  calculations. 
Olson,  Smith  and  Bauer  [54]  have  used  a  semi  empirical  method  to  obtain  appro¬ 
ximate  values  for  the  ^jf’s.  A  graph  was  set  up  similar  to  that  of  Bauer,  Fisher  and 
Gilmore  [55]  following  the  work  of  Hasted  and  Chong  [56] .  This  plot  included  all 
of  the  then  presently  available  versus  Wjf  values  for  one-electron  transfer  systems, 
98  points  from  both  experimental  and  theoretical  work.  It  was  found  that  by  plot¬ 
ting  reduced  quantities  the  data  scattered  about  common  curve  of  H*f  = 
versus  /? *  =  J  [(o  +  k)R^  ] .  Here  J  is  the  electron  affinity  of  B,  and  i  is  the 
effective  ionization  potential  of  A  in  eq.  (1).  When  parameterized  as 


//*f  =  /?;exp(-0.86f?p.  (26) 

eighty-three  percent  of  all  data  were  within  a  factor  of  3  of  this  curve.  The  range  of 
H*f  covered  was  from  10“*  to  10“'®  and  /?*  from  2  to  28.  Eq.  (26)  can  thus  be 
used  to  estimate  //^f. 

Another  method,  developed  by  Smirnov  [57]  can  be  used  for  the  calculation  of 
in  the  case  where  the  negative  iun  is  in  an  S  state,  it  is  much  more  rigorous  in 
nature  than  the  semi-empirical  scheme  of  Olson  et  al.  [54] .  and  employs  asympto¬ 
tic  expansions  for  the  wave  functions  of  the  negative  ion  and  the  excited  atom.  It 
is  represented  as 

=  70^  /1 2(  2«  2 )- 1  (4/e)l/7(a/?  ^)2/<>(  2/  +  1  )r-'  ( 1  /a  +  /  +  1 ) 

X  r->(l/o-0cxp[-(a  +  7)/?xl  ■  (27) 
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where  a  and  y  have  been  deHned  previously,  /  is  the  orbital  angular  momentum 
quantum  number  for  the  electron  in  the  excited  atom,  and/l^  =  2.65  for  H~ 

At  a  given  energy,  if  only  two  states  dominate  the  scattering,  eq.  (16)  may  be 
solved  in  closed  form  by  the  use  of  tabulated  integrals  [SI ,  S2j .  We  then  find  that 

Cif  =  47r«2(i +A£-/f)F3(X),  (28) 

where 

F3(X)  =  ;7e-^(l  -e-^’)2-^d2.  (29) 

in  which 

2u,  x/2n/?2//2(/?  )p'/2 

x= — 

F3(X)  is  effectively  an  integral  of  the  transition  probabilities  over  all  impact  param¬ 
eters.  The  important  range  of  crossing  distances,  iO  to  SO  a.u.,  corresponds  to  a 
range  of  AF  of  0.02  to  0. 1  a.u.  or  0.5  to  3  eV.  Thus  at  near  thermal  energies,  F  in 
eq.  (30)  becomes  negligible  compared  to  A£',and  X  approaches  a  constant  equal  to 
The  cross  section  is  then  domianted  by  a  v~-  dependence 

of  the  form 

=  AnRlAF  F^iX)IF  =  UnR^FjW/nv^  (31) 

where  F^CK)  is  constant.  This  v~  ^  tendency  is  evident  in  the  low  energy  experimen¬ 
tal  data  described  in  the  next  sections. 

For  energies  F  ^  AF,  Fj{\)  can  be  approximated  by  a  polynomial,  and  can 
be  parameterized  as 

Gjf  =/lV +fl/u  +  C  +  £>o.  (32) 

This  form  is  useful  for  extrapolating  merged  beam  data  to  lower  energies  for  the 
purpose  of  obtaining  thermal  energy  reaction  rates.  The  monoenergetic  reaction 
rate  is  a  =  v^Q,  and  thermal  rates  can  be  obtained  by  averaging  over  a  Boltzmann 
distribution. 

We  see  also  from  eq.  (3 1 )  that  at  low  relative  velocities  there  will  be  an  isotope 
effect  when  Q  is  plotted  against  o,.  At  high  velocities  Q  will  depend  only  on  for 
a  given  pair  of  reactants.  The  isotope  effect  will  be  shown  later  in  the  He'^  +  H~ 
and  He*  +  D  calculations. 

3.3.  Theoretical  difficulties 

Several  shortcomings  of  the  Landau-Zener  theory  and  the  curve  crossing  model 
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presented  above  should  be  recognized.  The  first  problem  is  associated  with  the  LZ 
theory  itself.  The  potentials  that  cross  are  assumed  to  be  linear  and  the  coupling 
matrix  element  is  assumed  to  be  constant  in  the  region  about  the  crossing  point 
R^.  These  approximations  can  lead  to  erroneous  results  if  the  potentials  deviate 
significantly  from  linear  forms.  There  are  also  problems  at  higher  energies  where 
the  Landau-Zener  theory  predicts  that  the  total  cross  sections  will  decrease  as 
f'*  -  instead  of  the  correct  £“*  dependence  |58,  59).  Recently  however.  Dubrov- 
skii  [60]  has  derived  a  correction  term  for  the  high  energy  LZ  dependence  which 
remedies  this  difficulty.  Also,  the  simple  LZ  theory  as  presented  above  does  not 
take  into  proper  account  the  possibility  of  tunneling  through  the  potential  bar¬ 
riers  to  allow  for  reaction  at  impact  parameters  whose  classical  turning  points  are 
at  internuclear  separations  greater  than  It  has  been  estimated  [42]  that  this 
effect  will  cause  the  calculated  lon-ion  neutralization  cross  sections  to  be  too  small 
h\  approximately  10^^. 

Most  likely,  a  good  sliare  of  the  difficulty  in  comparing  the  theoretical  results 
to  the  experimental  ones  does  not  simply  lie  in  the  inadequacies  of  the  LZ  theory 
but  in  other  aspects  of  the  problem.  One  source  of  difficulty  is  obtaining  coupling 
matrix  elements,  //jf.  that  are  more  accurate  than  a  factor  of  two  or  three.  This 
requires  detailed  ah  initio  potential  calculations  which  are  extremely  difficult, 
time  consuming,  and  costly  for  these  types  of  systems.  If  the  matrix  elements  are 
increased  by  a  factor  of  three,  there  can  for  many  cases  be  a  similar  increase  in  the 
calculated  total  cross  sections,  expecially  for  a  system  where  there  are  only  a  few 
product  states  available  and  therefore  few  curve-crossings. 

Another  possible  source  of  difficulty  in  previous  calculations  is  the  neglect  of 
rotational  coupling  between  reactant  and  product  states  of  different  symmetry, 
such  as  transitions  between  Z  and  II.  or  11  and  A  molecular  states.  For  the  IT^  +  H~ 
system,  it  has  been  found  [6 1  ]  that  the  inclusion  of  rotational  coupling  can  have 
an  effect  on  the  calculated  total  cross  sections  even  at  thermal  energies.  The  reason 
can  be  seen  by  examining  the  rotational  coupling  matrix  element  where  A is  the 
rotational  coupling  matrix  element  evaluated  at  R^.  If  we  look  at  the  impact  para¬ 
meter  that  corresponds  to  a  turning  point  at  R^.  cq.  (24)  and  substitute  it  into 
eq.  (3.^1  we  obtain  at  tiiermul  energies 

//,,(h  =  h,./?,)  =  uAJl  +(£^£)-'lV/?x  ■  (34) 

Thus,  if  is  of  die  order  of  0.1 .  even  at  thermal  energies  transitions  caused  by 
rotational  coupling  can  be  important  for  the  lighter  systems.  Noic  the  dependence 
in  eq.  (34)  on  the  reduced  mass  of  the  system. 

The  LZ  theorv  also  cannot  be  expected  to  apply  to  transitions  occuring  when 
there  is  no  curve  crossing.  Optical  measurements  performed  on  the  Na"*^  +  0~  system 
[43. 44]  show  an  appreciable  contribution  to  the  total  cross  section  from  a  product 
state  that  does  not  cross  the  Coulomb  potential.  This  measurement  indicates  that 
these  types  of  interactions  must  be  considered. 
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Even  with  ail  the  forementioned  difTiculties,  the  LZ  curve  crossing  theory  does 
display  an  energy  dependence  for  the  total  cross  sections  that  is  generally  in  good 
agreement  with  experimental  “merged-beam"  results.  Moreover,  the  predicted 
magnitude  of  the  cross  sections  are  within  a  factor  of  two  of  experimental  results 
for  most  species. 

3.3.  Close-coupled  calculations 

The  bulk  of  the  total  cross  section  calculations  for  atomic  species  of  reaction 
( 1 )  have  been  performed  using  the  L2  curve-crossing  model  because  of  its  ease 
of  application.  The  most  accurate  approach,  however,  would  be  to  perform  nu¬ 
merical  calculations  on  the  coupled  Schroedinger  equations  that  described  the 
nuclear  motion.  These  equations  must  be  solved  for  the  transition  probabilities 
of  each  possible  product  channel  at  every  orbital  angular  momentum  quantum 
number/.  The  probabilities  are  then  summed  to  obtain  the  cross  st:tion.  This  type 
of  calculation  requires  a  prohibi'ively  large  amount  of  computer  time  if  the  num¬ 
ber  of  possible  channels  exceeds  three  or  four.  As  a  result,  only  for  one  system. 

+  H~  has  a  calculation  of  this  type  been  attempted  [62] .  Here,  a/*  initio  po¬ 
tential  curves  and  coupling  matrix  elements  were  used  to  obtain  the  cross  sections. 
These  results  will  be  presented  in  the  next  section. 

One  approximate  dose-coupled  method  that  has  not  been  used,  but  appears  to 
be  ideally  suited  to  the  ion-ion  problem,  is  that  presented  by  Gordon  [63] .  In 
his  method,  the  potential  curves  are  approximated  by  linear  pieces  so  that  the  wave- 
functions  cart  be  represented  by  Airy  functions.  The  close-coupled  calculation  is 
then  performed  with  grids  equal  to  the  sizes  of  the  potential  pieces  instead  of 
being  determined  by  the  rapid  oscillation  frequency  of  the  wave  functions.  For 
the  ion-ion  case  where  the  product  potentials  are  almost  linear,  the  method  of 
Cordon  should  be  very  efficient  for  close<oupled  calculations. 

Classicul  close<oupled  methods  such  as  the  one  presented  by  Bates  and  Crothers 
[i)4]  could  also  be  applied  to  this  problem.  However,  the  difficulty  with  the  classi¬ 
cal  methods  is  that  an  “average”  trajectory  must  be  chosen.  In  the  case  where  there 
is  more  than  one  product  channel  this  becomes  extremely  difficult  for  low  energy 
collisions  (F  ^  100  eV).  At  higher  energies,  however,  the  classical  methods  may  be 
applied. 


4.  Atomic  ion  systems 
4.1.  If"  +H" 

A  relatively  simple  ion -ion  system  for  comparison  with  theoretical  calculations 
is  H'^  +  H~.  The  cross  section  for  this  reaction  has  been  measured  over  the  energy 
range  from  0.1 5  to  300  eV  at  SRI  [37] ,  and  from  125  to  5000  eV  by  Rundel.  Aitken 
and  Harrison  [46] .  The  results  are  shown  in  fig.  5.  The  two  measurements  are  in  ex¬ 
cellent  agreement  where  they  overlap,  and  in  combination  provide  a  measurement 
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t  ii;.  5.  Kxporimontal  and  Ihcoretiial  results  fur  ll'*'  +  H  .  •  •  Miiscloy,  Abcilh  and  Hcicrson  1 37 1 
*  -  Rundcl.  Althcn  and  llarriss>n  |46 1 ,  a  -  Bales  and  Lewis  |5 1 1 ,  b-  Dalyarno,  Vitlnr  and  Hlaiv 
cha/d  162|  .  c-Olson,  Peletson  and  Moseley  |4()| . 

Ilf  ihis  cross  section  spanning  an  energy  range  of  more  than  four  orders  of  magnitude. 

This  system  should  be  well  suited  for  a  test  of  theoretical  calculations  because 
there  are  only  two  electrons,  facili fating  afc  initio  calculations,  and  the  initial  and 
the  final  state  potentials  are  well  defined.  The  zero  order  potential  curves  for  this 
system  are  given  in  ftg.  (i.  A  number  of  l.Z  calculations  using  interaction  matrix 
elements  obtained  in  a  variety  of  ways,  have  been  made  on  this  system  and  the  re¬ 
sults  are  shown  in  fig.  .S.  The  earlie.st  calculation  by  Bates  and  Lewis  |.‘il  | ,  is  given 
by  curve  “a".  It  has  the  general  characteristics  of  the  measured  cross  section,  but  is 
lower  in  magnitude  by  about  a  factor  ol  three.  The  increase  in  the  cross  section  at 
low  velocities  is  caused  by  the  crossing  to  the  n  -  .3  stale  of  II.  and  the  maximum 
at  higlier  velocity  is  caused  by  the  n  =  2  crossing.  The  failure  of  the  calculation  to 
decrease  rapidly  enouglt  at  higli  velocity  is  expected  because  of  the  invalidity  of 
the  LZ  theory  at  these  high  energies.  At  higli  energies  the  LZ  theory  incorrectly 
predicts  that  the  cross  sections  should  decrease  as  I  'u  instead  of  \lir. 
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ATOMIC  UNITS  IiqI 
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INTERNUCLEAR  SEPARATION  R 
I'ii!.  6.  Zero  order  potenlial  curves  lor  H'*’  +  H  ^  •  II  +  II. 

Curve  "b"  represents  the  LZ  calculations  ot'Dalgarno,  Victor  and  Blanchard  |621 
which  were  based  on  more  accurate  potential  curves  for  the  H2  molecule.  This  re¬ 
sult  is  nearer  the  measured  cross  section  in  absolute  value  than  that  of  Bates  and 
Lewis  but  does  not  reproduce  the  velocity  dependence  satisfaclorily.  Recently. 
Browne  and  Victor  hare  extended  these  calculations  and  have  included  rotational 
coupling  between  initial  and  final  states.  They  found  |6I  ]  that  the  calculated  cross 
sections  arc  increased  substantially.  Apparently,  rotational  coupling  will  have  to  be 
included  in  any  detailed  calculation  on  an  ion-ion  neutralization  system. 

Curve  “c"  represents  the  LZ  calculation  of  Olson  et  al.  |40l  using  interaction 
nialri.x  elements  calculated  by  the  method  of  Smirnov  [.S7| .  This  calculalion  com¬ 
pares  quite  favorably  with  the  experimental  results  below  1 0  eV.  but  departs  from 
them  significantly  at  higher  energies.  The  Smirnov  formula  was  derived  using  asym¬ 
ptotic  wave  functkuis.  and  therefore  it  would  be  expected  th.it  interaction  matrix 
elements  for  crossings  at  large  distances  would  be  more  accurate  than  those  for 
smaller  distances.  Clearly  the  11  =  3  crossing  at  ^  36  a.u.  is  much  belter  repre¬ 
sented  than  the  11  =  2  crossing  at  ^  10  a.u. 

Recently.  Janev  and  Tancic  (651  have  calculated  the  11'*^  -t-  H  cross  section,  us- 
ign  the  technique  of  Dubrovskii  [55].  These  results  indicate  extremely  good  agree¬ 
ment  with  experiment.  However,  it  appears  there  is  a  numerical  error  in  the  partial 
cross  section  arising  from  the  11  =  2  state  so  that  the  higli  energy  cross  section  cal¬ 
culations  are  too  large  by  approximately  a  factor  of  3. 

It  thus  appears  the  LZ  calculations  are  able  to  reproduce  the  general  features 
and  approximate  magnitude  of  the  H*  +  H~  cross  section,  but  differ  from  the  ex- 
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perimental  results  in  a  number  of  significant  details.  Possible  reasons  for  these  dis¬ 
crepancies  will  be  discussed  later  in  this  section. 

Roy  and  Mukherjee  (66)  have  performed  a  close  coupled  calculation  on  H*  +  H 
over  the  lab  energy  range  from  0.5  to8  keV.  The  calculated  cross  section  is  in  agree¬ 
ment  with  experiment  at  8  keV,  and  is  about  309?  below  experiment  at  SOO  eV. 

This  represents  much  better  agreement  at  these  high  energies  than  any  of  the  LZ 
calculations. 

A  particularly  interesting  feature  of  the  experimental  +  H~  results  is  the  pre¬ 
sence  of  some  “fine  structure”  superimposed  on  the  broad  maximum  around  300 
eV.  In  addition  to  the  two  maxima  and  the  minimum  clearly  indicated  by  the  re¬ 
sults,  the  scatter  in  the  data  between  30  and  1 00  eV  indicates  that  there  could  be 
additional  structure.  Three  extrema  are  observed  that  are  rouglily  evenly  spaced  in 
reciprocal  velocity,  and  may  [32]  be  similar  in  origin  to  the  oscillations  seen  in  other 
cases  of  charge  transfer  or  energy  transfer.  Another  possible  source  of  this  structure 
is  the  variation  of  interaction  energies  with  internuclear  separation.  Bates  [58]  and 
Mordinov  and  Firsov  [67]  have  shown  that  the  inclusion  of  this  effect  in  an  LZ  cal¬ 
culation  can  possibly  lead  to  two  maxima  in  the  cross  section. 

4.2  N"" +0",0^ +0" 

The  experimental  data  on  N'*'  +  0~  of  Aberth  and  Peterson  [35] ,  grouped  ac¬ 
cording  to  energy,  avr  : 'd  and  converted  to  cross  sections,  are  presented  in  fig.  7. 
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Fig.  8.  E.xperiinentjl  |38|  and  (hcotclical  |40|  results  for  O*  +  o'~. 

along  with  LZ  calculation  of  Olson  et  al.  (40| .  This  calculation  is  substantially 
more  complicated  than  the  +  H~case  since  a  large  number  of  favorable  cros¬ 
sings  exist,  and  it  i.  necessary  to  include  ten  excited  N  states,  eight  excited  0  states, 
and  one  state  with  excitation  of  both  N  and  0.  The  calculations  are  further  com¬ 
plicated  by  the  fact  that  N*  +0  may  form  six  molecular  states,  requiring  calcula¬ 
tion  of  six  sets  of  cross  sections  weighted  by  the  appropriate  statistical  weight  fac¬ 
tor.  and  then  summed.  Details  of  the  calculation,  including  matrix  elements,  cros¬ 
sing  distances  and  energy  differences  for  all  states  considered,  are  given  in  ref.  [40] . 
It  can  be  seen  from  fig.  6  that  below  10  eV  the  theoretical  results  are  slightly  lower 
than  the  experimental  values,  but  the  energy  dependence  of  the  cross  sections  is 
well  reproduced. 

The  system  0"^  +  0”  has  also  been  studied  by  the  SRI  group  [38] .  The  results, 
shown  in  fig.  8  are  seen  to  be  very  similar  to  the  +  0~  resufis.  The  LZ  calcula¬ 
tions  [40]  are  somewhat  simpler  since  only  six  states  need  be  considered;  again 
the  theoretical  results  are  lower  than  the  experimental  values. 

4.3.  He^  +  H".  He^  +  D" 

Gaily  and  Harrison  [48]  have  investigated  He  +  H~  over  the  energy  range  from 
0.2  to  8  keV,  with  the  results  given  by  the  triangles  in  fig.  9.  The  SRI  group  has 
made  a  limited  number  of  measurements  on  He"*^  +  D~  [40] ,  represented  by  the 
circles.  At  relative  velocities  above  5  X  10^  cm/sec  the  two  cross  sections  should 


51 


J.T.  Moseley,  R.E.  Olson  and  J.R.  Peterson,  Ion  -ion  mutual  neutralization 


25 


BARYCENTRIC  ENERGY  (eV) 

Ht*  ♦  H- 

1  3  to  too  1000  10,000 

T — I - 1 — I - 1 — ^ - 1 - 1 - r 

♦  O' 


Fig.  9.  Expciiinental  results  on  He*  +  H”  148|  (A)  and  He*  +  D”  |40|  (o),  and  theoretical  re¬ 
sults  for  both  systems  |40|  using  scmi-empiiical  matrix  elements  (dashed  lines)  and  the  Smir¬ 
nov  formula  |54|  (solid  lines). 

be  equal.  The  agreement,  as  evidenced  by  fig.  9,  is  satisfactory.  Measurement  at 
lower  energies  were  not  made  because  at  the  time  these  measurements  were  made 
the  SRI  merged  beam  apparatus  was  unable  to  merge  ions  of  mass  ratio  2  at  very 
low  relative  energies.  This  limitation  has  since  been  removed  (38]. 

The  LZ  calculations  of  Olson  et  al.  [40]  are  also  shown  in  fig.  9.  Here,  the  dash¬ 
ed  lines  correspond  to  the  matrix  elements  calculated  by  the  Smirnov  formula  and 
the  solid  lines  refer  to  matrix  elements  estimated  by  the  semiempirical  formula 
[54] .  The  discrepancy  in  energy  dependence  above  5  keVean  again  be  explained 
by  the  invalidity  of  the  LZ  theory  at  these  high  energies.  At  low  energies  the  cross 
section  calculations  predict  the  isotope  (mass)  effect  mentioned  in  section  3. 

The  structure  in  the  He*  +  H~  measurements  is  actually  more  pronounced  than 
it  appears  in  fig.  9.  This  structure  occurs  in  the  same  velocity  range  as  the  structure 
in  the  H*  +  H~  cross  sections,  and  most  likely  has  the  same  origin. 

4.4.  Na*  +  0~ 

The  reaction  Na*  +0  -♦  Na  +  0  was  first  studied  by  Weiner,  Peatman  and  Berry 
[43,44] ,  who  observed  photons  emitted  by  the  Na  products,  as  described  above  in 
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Fig.  10.  Zero  order  poicniial  curves  for  Na^  +  0~  -*  Na  +  O. 

section  2.2.  The  zero  order  potential  curves  for  this  system  are  presented  in  fig.  10. 
Measurements  were  made  on  photons  from  the  Na  3^0  3^?,  3-S.  and  4^P  -►  3^S 

radiative  transitions.  The  cross  sections  obtained  from  the  3^0  -♦  3^P  transitions 
are  shown  in  fig.  1 1 .  They  are  the  largest  yet  reported  for  a  mutual  neutralization 
reaction.  These  results  were  suprising  because  the  Na^  +  0”  Coulomb  potential 
crosses  the  Na(3^D)  +  0(^P)  curve  at  /?  >  200  A  and  the  LZ  theory,  using  the  semi- 
empirical  matrix  elements  [54] ,  predicts  a  negligibly  small  cross  section  for  a  curve 
crossing  at  such  a  large  internuclear  separation.  Even  more  surprising  from  the  curve¬ 
crossing  viewpoint  is  the  observation  of  a  cross  section  of  about  100  A  for  '.he  pro¬ 
duction  of  Na(42p),  for  which  the  reaction  is  0.1  eV  endothermic.  Such  jn  endo¬ 
thermic  reaction  could  depend  on  a  curve  crossing  with  the  repulsive  portion  of 
the  Coulomb  potential  at  ^  3  A,  but  this  cross  section  would  have  an  upper  limit 
of  less  than  nR^,  or  about  30  A^. 

The  cross  section  obtained  by  observing  the  3P-v  3S  transition  should  be  greater 
than  that  for  the  3D  -►  3P,  since  it  will  include  all  the  3D  -♦  3P  cascade  contributions 
plus  any  direct  contribution  from  the  Na(3P)  +  0(^P)  crossing  at  ~  17ao'  In  addi¬ 
tion.  it  should  be  less  than  the  total  neutralization  cross  section  since  it  will  not  in- 
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Fig.  1 1.  Experimental  results  ublaincd  from  observation  of  the  Na  3^D  -»  3*P  transition  |44|. 

dude  contributions  from  the  observed  4P-*-  3S  transition,  or  from  the  Na(3S)  + 
0(*D)  crossing  at  ~  Ibog. 

The  total  cross  section  has  been  rr.easured  by  the  SRI  group  (38| ,  and  is  shown 
on  fig.  1 2,  along  with  the  3P  -♦  3S  results  of  Weiner  et  al.  and  a  Landau-Zener  cal¬ 
culation  performed  by  Olson  (68| .  The  SRI  total  cross  section  measurements  are 
believed  to  be  accurate  to  within  ±  30%,  and  the  absolute  uncertainty  in  the  3P-3S 
measurements  is  placed  at  a  factor  of  2.  The  3P  -►  3S  cross  section  are  about  three 
times  larger  than  the  total  cross  sections;  thus  even  considering  these  error  limits 
a  discrepancy  exists  between  the  two  measurements. 

The  magnitude  of  the  3P  -»  3S  cross  section  is  also  inconsistent  with  that  of  the 
3D  -►  3P  cross  section.  Rather  than  being  larger,  it  is  in  general  smaller,  and  at  5  eV 
is  about  a  factor  of  three  smaller. 

As  we  have  discussed  earlier,  determination  of  absolute  cross  sections  by  this  op¬ 
tical  technique  is  quite  difficult,  and  perhaps  this  inconsistency  arises  from  a  cali¬ 
bration  error.  In  any  case,  it  is  our  belief  that  the  total  cross  section  measurements 
of  the  SRI  group  give  a  more  reliable  result  for  the  absolute  magnitude  of  these  cross 
sections.  Nonetheless,  the  photon  measurements  of  Weiner  et  al.  are  of  considerable 
importance,  since  even  if  the  actual  3D  -►  3P  cross  section  is  an  order  of  magnitude 
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I'ig.  1 2.  Results  for  the  totul  cross  section  Na'^  +  O  -  Na  +  O  1 38|  (a),  from  observation  of 
the  Na  3^P  3^S  transition  |44|  (*),  and  from  a  LZ  calculation  |68)  (dot  -dash  line), 


smaller  than  reported,  the  validity  of  the  simple  LZ  theory  to  predict  product  states 
is  made  questionable.  Further  investigations  of  emitted  photons  should  shed  con¬ 
siderable  light  on  the  mechanism  of  ion-ion  interactions. 

The  results  of  the  Na^  +  0~  measurements  encourages  a  more  detailed  considera¬ 
tion  of  the  final  states  of  the  reactioi  'han  was  presented  in  the  discussion  of  the 
simple  LZ  theory  in  section  3.  That  discussion  included  only  states  that  interact  on 
the  descending  part  of  the  Coulomb  curve.  Collision  channels  that  remain  on  the 
Coulomb  curve  and  rcuch  tlie  repulsive  wall  of  the  potential  at  ^  3  A  can  also  be 
considered.  Interactions  of  the  Na(42p)  +  Of^P)  states  are  then  not  only  possible 
but  should  be  expected.  Oscillations  such  as  those  observed  on  the  cross  sections 
for  Na(3^D)  and  Na(3^P).  which  are  roughly  spaced  as  1/uand  are  out  of  phase 
with  one  another,  can  result  from  transitions  from  the  initial  Colomb  state  to  in¬ 
termediate  states  that  are  mixed  at  larger  separations  at  an  avoided  curvfe  crossing 
via  the  Rosenthal  mechanism  [69].  However, calculations  based  on  this  model  [68] 
of  reaction  on  the  repulsive  wall  of  the  interaction  potentials  yield  cross  sections 
about  an  order  of  magnitude  smaller  than  those  reported  by  Weiner  et  al. 


55 


J.  T.  Moseley,  R.l\.  Olson  and  J.K.  Peterson,  Ion  ion  nnitiial  neiilralizatUin  29 


lablc  1 

Roaclion  lalc  aK-HKicuts  al  300  K  aloinjc  ions  a( 300  K)  10'^  ciii  ’  m-i.- 


System 

lApcrimenlal 

11^  +  ir 

3.9  •  2.1 

1,5-' 

4.tl'’ 

1.2  |51| 

+  o' 

2.6  .*  (1.8 

1.8-' 

+  0 

2.7  1  1,3 

l.l-' 

'(1.8  |50| 

Na'^  +  0 

2.1  1  Id 

n.7'^ 

iio^  +  ir 

5.7'' 

7.3'' 

11C++ 

_ 

4.7^ 

5.7'’ 

^  Total  cross  sections  calculated  usiiiu  s’O.  (26). 

L  '  * 

“  Total  cross  sections  calculated  usinjtcq.  (27). 

'•  It  sliould  be  noted  tliat  the  cloclroncirativivty  ofO  used  in  this  ref.  150|  seas  2.2  cV,  whereas 
the  value  now  accepted  is  1 .47  eV.  Tse  of  this  value  inipht  suhstaiitially  al'Icet  the  calculated 
rate  constant. 


4.5.  Thermal  energy  reacthni  rate eoej'ju  icnis 

A  primary  interest  iti  ilic  iim  -ion  neoirali/alion  reaction  is  its  imponatice  in 
iotii/ed  gases,  Ftir  such  applications  the  parameter  orimpoilance  is  ilie  reaclioti 
rate  coelTicicnt  a  at  thermal  temperatures.  This  rate  can  be  estimated  from  the 
SRI  data  in  the  tnanner  discussed  below. 

As  discussed  in  section  3,  the  LZ  fortuula  can  be  approxitnaled  at  low  energies 
by  the  expansion 

Q=A/vj  -i-B/u^+C  +  Dv^.  (32) 

Experimental  data  can  be  easily  fit  to  such  parameterization,  and  extrapolated  to 
lower  energies.  The  low-energy  temperature  depend ;ncc  of  the  ihermal  rate  coef¬ 
ficient  Q  =  <(3Uj>  can  then  be  calculated  [381  by  Boltzmann  averaging  v^Q.  Error 
limits  on  o  can  be  obtained  from  the  uncertainty  in  the  fitted  curve  and  in  the  data. 
The  dashed  curves  through  the  data  in  figs.  7  and  8  represent  such  parametcri/a- 
tions. 

Table  1  gives  reac  tion  rales  at  300°K  for  +  H  .  hT  +0  .0^+0  ,  and 
Na^  +  0~  as  determined  from  the  data  presented  here.  In  addition,  theoretical  LZ 
values  of  a  arc  given  for  these  systems  and  for  Hc'^  +  and  fie'^  +  D". 

4.6.  Discussion 

From  the  measurements  and  calculations  presented  here,  a  number  of  generali¬ 
zations  can  be  made  about  atomic  ion  -ion  mutual  neutralization.  Reaction  rates  at 
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llK'rm;il  energies  for  siieli  systems  can  be  expected  to  exceed  10"^  cm^/sec.  The 
cross  section  at  energies  below  a  few  tenths  of  an  eV  can  be  expected  to  vary  es¬ 
sentially  as  1//:',  Cross  sections  calculated  using  the  LZ  theory  are  usually  smaller 
than  are  observed  experimentally,  but  are  within  a  factor  of  2  or  3  of  the  measur¬ 
ed  cross  sections,  and  they  predict  the  general  velocity  dependence  of  the  cross 
sections  at  low  to  intermediate  velwities. 

The  general  tendency  of  the  LZ  calculations  to  be  somewhat  smaller  than  ex¬ 
perimental  results  could  arise  from  a  number  of  factors.  The  ions  produced  in  the 
laboratory  might  not  all  be  in  their  ground  states,  and  hence  a  coupling  trr  more 
states  could  result,  giving  rise  to  a  large  cross  section.  The  LZ  theory  may  not  pro¬ 
perly  account  for  the  range  of  y?  over  which  an  interaction  can  occur.  Further,  the 
LZ  theory  clearly  dues  not  allow  significant  interaction  at  large  or  small  distances, 
but  both  regions  may  be  important.  Evidence  for  this  is  found  in  the  observation 
of  Weiner  et  al.  of  large  cross  sections  for  the  Na^  +  0  (3D  3P)  transition,  even 

thougli  the  crossing  is  at  nearly  300  A.  and  for  the  4P  -*■  3S  transition,  even  though 
no  crossing  occurs  with  the  Na(4P)  state  until  the  repulsive  part  of  the  potential 
is  reached.  Similar  interactions  might  account  for  most  of  the  discrepancies  between 
theory  and  experiment  noted  in  this  section.  For  example,  the  failure  at  low  ener¬ 
gies  of  the  most  rigorous  calculations  of  H'*’  +  H~  (curve  b  of  fig.  5)  may  be  due  to 
the  failure  of  the  LZ  formalism  to  properly  account  for  the  crossing  at  large  dis¬ 
tances  to  the  I.  =  4  levels  of  hydrogen.  More  theoretical  work  is  needed  to  under¬ 
stand  the  importance  of  transitions  between  curves  which  are  nearly  degenerate 
over  long  distances  and  to  investigate  the  importance  of  close  encounters.  Addi¬ 
tional  experimental  measurements  of  the  final  states  of  the  resulting  neutrals  are 
also  needed  to  verify  the  theoretical  work. 


5.  Cross  section  calculations  for  molecular  systems 

It  becomes  very  difficult  to  apply  the  LZ  method  presented  in  section  3  to  the 
system 

A^  +  B  -*A  +  B  +  At'  ( 1 ) 

when  A  or  B  or  both  are  molecules.  Now.  the  number  of  curve-crossings  with  the 
reactant  Coulomb  potential  can  easily  reach  into  the  hundreds  due  to  the  large  num¬ 
ber  of  excited  electronic  states  and  vibrational  levels  available.  Even  if  the  coupling 
matrix  elements  could  be  evaluated  at  each  curve  crossing,  the  computer  time  need¬ 
ed  to  solve  for  the  transition  probabilities  via  the  simp!;  LZ  method  becomes  pro¬ 
hibitive.  Therefore,  some  other  approach  must  be  used. 

One  reasonable  assumption  is  that  for  the  more  complicated  systems  the  physi¬ 
cal  situation  may  be  approximated  using  a  high  density  of  crossing  states.  Then  by 
using  the  Landau-Zener  method  and  semi-empirical  coupling  matrix  elements,  a 
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critical  crossing  distance  can  be  calculated  by  assuming  an  absorbing-sphere  model, 
i.e.,  unit  probability  exists  for  reaction  within  This  approach  is  discussed  in  detail 
in  ref.  [42] ,  and  will  be  summarized  here.  To  determine  the  critical  distance 
within  which  there  is  reaction,  consider  the  application  of  the  two-channel  Landau- 
Zener  method  discussed  in  section  3.  In  this  approximation  the  total  cross  section 
is  given  by  eq.  (28)  where  for  the  special  case  of  thermal  energy,  X,  from  eq.  (30), 
approaches  a  constant  given  by  \/2itR^^ two-channel  case, 
the  integral  of  the  transition  probabilities,  f3(X),  eq.  (29),  has  a  maximum  value 
of  0. 1 1 3  when  X  =  0.424. 

If  we  now  extend  these  results  to  the  multichannel  case  where  a  large  number 
of  final  states  are  available,  we  may  expect  reactions  for  all  trajectories  that  ap¬ 
proach  within  some  critical  internuclear  separation  R^.  Using  the  two-channel  case 
as  a  guide,  R^.  will  be  slightly  larger  than  the  R  value  that  satisfies  X  =  0.424.  The 
total  cross  section  from  eq.  (28)  will  then  be 

G  =  7r/?2|l+(«^,£^“'|,  (35) 

To  determine  X,  we  must  know  its  dependence  on  R.  The  coupling  matrix  ele¬ 
ment  //[  2  is  the  only  unknown,  and  from  previous  work  on  one-electron  transfer, 
we  have  been  able  to  parameterize  it  in  terms  of  the  effective  ionization  potential 
of  the  reactants  and  products.  The  relation  was  given  in  eq.  (26)  of  section  3.  This 
functional  form  forW|2(/?x)  is  valid  only  at  large  distances,/?  ^  lOflo,  where  the 
electron  is  transferred  through  the  potential  barrier  determined  by  the  exponen¬ 
tial  tails  of  the  wave  functions  for  the  reactant  and  product  states  and  for  values 
of  /j  and  If  greater  than  ~  0.4  eV. 

By  knowing  the  electron  detachment  energy  of  the  negative  ion  and  assuming 
that  there  is  a  high  density  of  final  states  available  for  reactions,  we  may  determine 
7  as  a  function  of  R  using  the  relationship 

/f  =  /(+/?-'.  (.36) 


Using  eq.  (26),  it  is  an  easy  matter  to  obtain  X,  eq.  (30),  as  a  function  of  R.  From 
numerical  calculations  that  include  large  numbers  of  final  stater,  we  find  that  X  ap¬ 
proaches  0. 1 5  as  the  density  of  the  states  approaches  a  continuum.  The  value  X  = 
0.424,  which  is  correct  for  the  two-channel  case,  tends  to  underestimate  the  cross 
section  by  about  10^.  The  use  of  X  =  0.1 5  to  determine  R^  has  also  been  verified 
by  Landau-Zener  calculations  employing  a  large  number  of  product  channels. 

Because  of  the  Coulomb  attraction,  the  product  EQ  approaches  a  constant  at 
very  low  (i.e.  thermal)  energies  in  eq.  (35),  so  that  the  reaction  rate  a,  which  is  an 
average  over  a  Maxwellian  distribution  of  ion  velocities,  becomes 


u^du  2(u^Q) 


(37) 
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For  molecular  systems,  the  coupling  matrix  elements,  eq.  (26),  should  be  mod¬ 
ified  by  the  Franck-Condon  factors  4  that  represent  the  overlap  of  the  vibrational 
energy  levels,  to  read  (50] 

//*2  =q^'^R*  exp(-0.86«*) .  (38) 

The  Franck-Coiidon  factors  are  quantities  less  than  unity  so  not  including  them  in 
the  calculation  should  provide  an  upper  limit  on  the  cross  section.  However,  for 
the  cases  studied,  even  if  all  were  set  equal  to  0.1,  the  cross  sections  would  be 
decreased  by  only  about  20%..  Therefore  the  results  are  not  strongly  affected  by 
ignoring  the  Franck-Condon  factors. 

The  simplified  treatment  above  considers  cases  when  a  high  density  of  final 
states  exists  in  the  region  of  most  favorable  curve  crossings.  Cross  sections  from 
these  formulae  thus  represent  upper  limits  to  the  actual  values.  On  the  other  hand 
cross  sections  calculated  in  this  manner  are  probably  not  much  larger  than  the  ac¬ 
tual  values  even  when  the  density  of  final  states  is  low,  since  each  reaction  channel 
can  then  have  a  relatively  high  probability  due  to  the  lack  of  competition  with 
neighboring  channels.  When  the  number  of  reactant  states  is  larje  (n  ^  1 0),  the 
cross  section  is  not  strongly  dependent  on  the  number  of  states.  The  cross  section 
is  then  found  to  approach  closely  the  absorbing  sphere  value. 

In  the  absorbing-sphere  model  described,  the  cross  section  is  dependent  only 
on  the  value  of  the  electron  detachment  energy  of  the  negative  ion  and  the  reduced 
mass  of  the  system.  A  similar  result  is  obtained  in  the  work  of  Radtsig  and  Smir¬ 
nov  [70]  who  use  a  different  approach  to  this  problem. 

In  a  given  system,  as  the  electron  detachment  energy  is  decreased  because  of 
vibrational  or  electronic  excitation,  the  reaction  rate  increases.  Physically,  this  cor¬ 
responds  to  the  outer  electron  011  the  negative  ion  being  more  loosely  bound,  mak¬ 
ing  it  possible  for  transfer  to  occur  at  larger  intermolecular  separations.  If  the  negative 
ion  is  in  an  excited  electronic  or  vibrational  state,  we  would  predict  that  the  reaction 
rate  will  generally  be  gre.iter  than  when  the  negative  ion  is  in  its  ground  state.  This 
problem  usually  arises  with  molecular  negative  ions,  but  may  also  occur  with  an 
atomic  negative  ion  system  if  the  negative  ion  has  a  metastable  electronic  state. 

In  order  to  account  for  transitions  that  occur  in  a  range  around  RJ^,  the  Landau- 
Zener  transition  probabilities  can  be  compared  with  the  close-coupled  values  tabu¬ 
lated  in  reduced  form  by  Delos  [71  j .  It  is  found  that  the  cross  sections  obtained 
by  the  Landau-Zener  method  are  below  the  close-coupled  results  by  approximate¬ 
ly  10%.  This  factor  is  included  in  the  absorbing  sphere  calculations  to  be  presented 
here. 

The  energy  dependence  of  the  cross  sections  may  be  calculated  using  eq.  (30) 
for  the  evaluation  of  X.  In  fig.  13  are  shown  the  cross  section  curves  predicted  by 
the  absorbing  sphere  model  for  the  negative  ions  studied  here.  The  curves  for  each 
system  are  essentially  independent  (within  10%)  of  the  positive  ion  mass  when  it  is 
within  a  factor  2  of  that  of  the  negative  ion. 
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As  expected,  the  system  with  the  lowest  electron  affinity,  O2,  has  the  largest 
cross  section  at  a  given  energy.  Likewise,  the  system  with  the  largest  electron  affin¬ 
ity,  NO3,  has  the  smallest  cross  section.  If  the  negative  ion  is  in  an  excited  level 
above  its  ground  state,  we  would  expect  the  observed  cros  sections  to  be  larger 
than  those  predicted  here. 

A  few  comments  should  be  made  concerning  the  significance  of  the  dissociation 
reaction 

A-fBC^-A  +  B+C-*- A£  (39) 

since  this  may  be  the  source  of  some  of  the  theoretical  and  experimental  differences 
in  the  molecular  systems.  Bates  and  Boyd  (52)  have  presented  Landau -Zener  for¬ 
mulas  for  this  process,  and  we  find  that  the  absorbing  sphere  model  will  be  equally 
applicable  to  this  type  of  reaction.  The  only  difference  between  the  theory  for 
reaction  ( 1 )  and  that  for  reaction  (39)  is  in  the  energy  dependence  of  the  cross  sec¬ 
tions.  At  low  energies,  in  the  region  of  practical  interest  for  reaction  rates,  both 
reactions  ( 1 )  and  ( 39)  will  have  a  cross  section  that  varies  approximately  as  ' . 

At  higher  energies,  the  cross  sections  for  reaction  ( 1 )  continue  to  decrease,  but  at 
a  rate  slower  than  E  ’ ,  until  energies  above  about  I  keV  are  reached,  when  the  de¬ 
pendence  again  approaches  f .  The  cross  sections  for  reaction  (39)  at  higli  ener¬ 
gies.  however,  approach  a  constant  value. 

The  cross  section  behavior  at  high  energies  therefore  will  depend  on  whether  or 
not  there  are  dissociating  product  channels  available  for  reaction.  We  may  predict 


Fiji.  13.  Ion -ion  mutual  neutralisation  cross  sections  taleula  ted  using  the  absorbing  sphere  mod¬ 
el  (42| .  In  this  model  the  cross  section  has  only  a  very  weak  dependence  on  the  mass  of  the  posi¬ 
tive  ion. 
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th;it  reuL'tiuii  (3‘))  will  be  inipiirtant  when  neutral  slates  iifthe  positive  nioleeular 
ion  have  curve  crossings  in  the  iinpoiiani  10  50U|,  region,  which  lead  to  dis¬ 
sociation.  These  conditiotis  will  occur  for  dissociated  neutral  states  that  lie  below 
the  ground  ion  slate  by  about  0..S  to  T.O  eV.  plus  the  electron  detachment  energy 
of  the  negative  ion.  Sii  far.  the  energy  dependence  of  the  experimental  data  on  the 
niolecniar  ions  tends  to  indicate  that  the  dissociation  reaction  is  not  a  dominant 
process  for  those  systems  measured. 

With  regard  to  the  vibrational  effects  in  molecular  ions,  the  most  I'avoied  elec¬ 
tronic  transitions  (which  occur  at  large  distances)  will  be  to  states  with  the  greatest 
Franck  Condon  I'a  .tors.  The  extent  to  which  vibr;  lional  energy  is  changed  during 
the  collision  thus  depends  largely  on  the  sliapes  and  portions  of  the  potential  energy 
curves  of  the  incident  ions  and  product  neutrals.  It  is  reasonable  to  expect  that 
some  vibrational  excitation  is  present  in  most  bound  product  molecules. 

Probably  the  most  signiticani  factor  to  note  conceining  the  absorbing  sphere 
model  Is  the  large  reaction  rate  dependence  on  the  electron  deiachmeni  energy . 

If  the  electrmt  detachment  energy  is  decreased  because  of  populated  excited  slates 
of  the  negative  ion.  the  cross  sections  increase.  If  the  excitation  is  to  the  upper 
levels,  this  increase  becomes  signiticani.  For  upper  atmosphere  chemistry .  this  fad 
makes  it  desirable  to  know  the  excited  state  populations  of  the  reactants  foi  the 
ion  ion  recombination  reaction. 


6.  Simple  molecular  systems 

6./.  m;  +  d“ 

The  11,  +1)  reaction  represents  a  simple  molecular  ion  system  and  consequently 
provides  an  impiMtanl  test  for  any  theory  describing  neulrali/alion  reaction  involv¬ 
ing  molecular  ions.  The  results  of  measurements  on  this  system  1 72|  are  given  in 
fig.  14.  The  solid  curve  represeiils  an  absorbing  spheie  calculation  |fi8| .  There  is  ap¬ 
parently  some  oscillatory  structure  in  the  experimental  data  above  5  eV.  alihouglt 
the  magnitude  ol  the  errors  casts  some  doubt  on  this  conclusion.  The  absorbing 
sphere  calculation  of  the  rale  coefficieirl  at  .300°  K  yields  the  value  of  ( 8.5 1  1.5).* 
Fxlrapolation  of  the  experimenial  results  y  ields(4.7  +  1.5). 

ft.  2.  NO*  +0".  0,  +0“ 

Results  for  the  reactions  NO*  +  0~  and  O*  +  O  1.38)  arc  given  in  figs.  I  5  and 
lb.  The  data  here  are  plotted  as  the  product  of  the  relative  velocity  u,  and  the  cross 
section  Q.  This  product  is  effectively  a  mono-energetic  reaction  i.ile.  The  energy  de¬ 
pendence  of  these  rales  are  typical  for  ion -ion  neuirali/ation  involving  complex 
ions.  The  solid  curves  in  the  figures  represent  a  least-squares  lit  of  the  data  to  the 
functional  form  ofcq.(.3Z). 

*  All  rate  eoerfieiciils  w  ill  be  expressed  in  lire  iiniis  in"’  irn'*/see.  and  these  iiniis  will  be  omitted 
in  llie  lexi. 
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BARYCENTRIC  ENERGY  (tV) 


Mg.  14,  l-xpaimcntal  |72|  anil  thcorclital  |58|  usullsim  II2  +  D 


BARYCENTRIC  ENERGY  |tV| 
ai  0.S  1  3  6  10  16  20 


I  ig  15  Kxporimcnlal  roulix  (3X|  on  NO*  +  O  . 

It  is  interesting  to  compare  these  reactions  with  the  atomic  ion  reactions  that 
have  been  investigated  lor  0  ,  Fig.  1 7  shows  the  cross  sections  for  all  live  prisitive 
ions  that  have  been  studied  with  0~.  All  of  the  reactions  show  a  very  similar  en¬ 
ergy  dependence,  and.  except  for  NO*,  have  very  nearly  the  same  magnitude. 
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BARYCENTRtC  ENERGY -  iV 
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l  ig.  16. 1’xpciimcnliil  results  1 38)  on  +  O 


l  ij:.  1  7.  Results  lor  the  Meutr.ili/jtion  of  N*.  O*.  Na*.  ot.  jnU  NO*  b\  O’.  The  .sohJ  eurve  is 
lioni  an  absorbini;  sphere  ealeiilatton  |421 


The  solid  curve  represenis  an  absorbing  sphere  caleulation  by  Olson  |42| .  As  has 
been  discussed,  for  this  model  the  cross  section  is  dependent  only  on  the  value  of 
the  eleclion  detachment  energy  of  the  negative  ion  and  the  reduced  mass  of  the  sys- 
icin.  It  is  striking  that  such  a  simple  model  yields  quite  good  agreement  with  the 
experimeiilal  results,  and  that  in  lad  the  ross  sections,  except  for  NO*,  are  very 
nearly  equal. 

The  experimental  cross  sections  were  extrapolated  to  lower  energies,  and  the  reac¬ 
tion  rates  at  300°  k.  were  calculated  in  the  manner  described  earlier.  The  results,  for 
NO*  and  O,  respectively,  were  (4,4  ±  2.0)  and  (1,0+  0.4),  Absorbing  sphere  calcula¬ 
tions  predicl  ( 1 .4  +  0.6)  for  both  of  these  reactions. 
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6..?.  I,  +  I 

Recunibinution  lor  iodine  ions  formed  in  a  gas  has  been  studied  at  temperatures 
near  300°  K.  Yeung  1 1  ]  observed  the  change  in  tiie  dielectric  constant  in  iodine 
vapor  due  to  the  decay  of  ions  following  a  discharge.  Greaves  125|  used  the  same 
technique,  but  added  a  mass  spectrometer  to  identify  the  ions  present.  He  also  coi- 
rected  a  calibration  error  in  Yeung's  e.xperiinent  which  had  caused  Yeung's  report¬ 
ed  values  to  be  too  small  by  a  factor  of  10.  Carlton  and  Mahan  [281  formed  ions 
by  photo-ionization  and  used  charge  collection  to  monitor  the  decay  of  the  ion  con¬ 
centration.  The  values  obtaitied  for  a  at  room  temperature  by  these  three  e.xperi- 
ments  are,  respectively .  ( 1 .47  ±  0.07),  ( 1 .22  ±  0.03)  and  ( 1 .4.S  +.  0.5).  .Although 
the  first  two  values  do  not  agree  within  the  quoted  errors,  the  overall  agreement 
for  these  three  expeiiments  is  quite  good.  However,  it  has  been  pointed  out  |2I  ] 
that  the  condition  that  diffusion  effects  be  negligible  is  only  marginally  satisi'ied 
by  all  of  these  expeiiments.  and  as  a  result,  each  of  these  results  may  be  high  by  as 
much  as  20  to  30'/'. 

Yeung  and  Greaves  investigated  the  temperature  dependenee  of  ihe  Ij  *  l~  reac¬ 
tion.  The  recombination  rate  was  observed  to  decrease  with  increasing  temperature, 
as  expected.  .Although  the  temperature  range  was  quite  limited,  both  results  appear 
to  indicate  a  decrease  in  die  cross  section  substantially  fastci  than  the  expected 
|f/;'(or  l/D  dependence. 


7.  Complex  molecular  systems 

Recombination  systems  in  which  both  positive  and  negative  ions  are  molecular 
are  quite  complex  since  a  large  number  of  potential  surfaces  will  be  involved  in  the 


RELATIVE  ENERGY  eV 


I  iv  18  Result'  I'T  l!u'  aein r.iii/alion  iil  ot .  NO  .  .iiut  Ni  bv  0; .  1  lie  solid  eiirve  is  Ironi  an  ab- 
surbinc  spliere  s.ileiil.iii  'ii  IbS] 
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reaction.  Very  little  is  known  about  most  of  these  poteniials.  and  a  wide  variety  ol' 
final  states,  both  dissociating  and  non-dissociating,  are  likely  to  be  involved.  It  is 
ions  ot  this  type,  however,  that  are  responsible  lor  the  behavior  of  most  ioni/ed 
gases  ot  juactical  importance.  We  shall  discuss  reactions  that  are  grouped  accord¬ 
ing  to  the  negative  incident  ion. 

7.1.  with  (>7 

Measurements  have  been  made  on  the  neutrali/.ation  of  07  with  |.'5| .  Ot, 
l-^)) .  and  NO^  (7.3| ,  The  re.sultsarc  .slrown  in  fig.  IK.  The  solid  curve'is  from  an 
absorbing  .sphere  calculation  |6K| .  The  NO^  and  cross  sectiotis  are  ver}'  tiearlv 
equal,  whtle  the  cross  section  is  smaller,  particularly  at  lower  energies. 

h.Ntrapolatioti  ot  the  higher  energy  +  O,  results  showti  in  fig.  IK  yields 
(4.2  ±  1.3)  .iitd  an  absorbitig  sphere  calculation  yields  ( 2.4  ±  O.K).  A  tneasurement 
of  the  -I-  0,  neittrali/ation  rate  at  thermal  energy  has  been  reported  by  Hirsh 
(ttid  histier  1 74| .  They  obtaitted  a  value  of  ( 1 .0  ±  0. 1 ).  Possible  reasotts  for  those 
differettces  will  be  disettssed  later  in  this  section. 

7.2.  .\'ciiiriili:atit)n  with  N07 

Neitirali/atioti  rates  have  been  measured  for  +  N07  atid  NO^  -r  N07  |.V)| 
over  the  eitergy  range  frotn  0.1  5  to  200  eV  and  are  shown  in  figs.  1‘)  atid  20.  The 
rates  and  etiergy  dependettcies  of  these  two  reactions  are  qitite  differetit  at  ettei- 
gies  above  5  eV.  but  below  this  energy  the  rates  are  nearly  equal,  h.xtrapolation 
of  these  results  yields  a  rate  at  .^OO”  K  of  (4.1  ±  1 .3)  for  and  ( 5. 1  ±  K.'^lfor 
NO^.  Very  recent  measurements  at  SRI.  as  yet  unpublished,  yield  a  300°  K  rate 
for  N'li  f  NO,  of  ( 1 .3  i  0..S). 
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I'ig.  20,  Kxpciimcntal  ivsults  |39|  on  NO^  +  NO^. 

A  thermal  energy  (300°  K)  rate  of  ( 1 .75  ±  0.6)  has  been  measured  for  NO^  + 
NO2  by  Eisner  and  Hirsh  (26| .  Mahan  and  Person  (27)  have  reported  a  rate  of 
(2.1  ±  0.6)  for  ions  believed  to  be  primarily  NO''^  and  NOj.  Absorbing  sphere  cal¬ 
culations  yield  ( 1 .2  ±  0.3)  for  both  O*  and  NO*,  and  ( 1 .3  ±  0.3)  for  N^. 

7.3.  Neutralization  with  NOJ 

Eisner  and  Hirsli  [26)  have  reported  a  thermal  rate  of  (0.34  ±0.12)  fur  NO*  + 
NOj.  Tins  is  the  smallest  value  for  a  thermal  energy  rale  that  has  been  reported.  Ab¬ 
sorbing  sphere  calculations  yield  a  value  of  ( 1 . 1  ±  0.3).  and  merged  beam  measure¬ 
ments  [72]  yield  a  value  of  3  at  0.1 5  eV,  and  extrapolate  to  (8.1  ±  2.3)  at  thermal 
energy.  Clearly  there  is  a  substantial  discrepancy  in  these  values. 

The  merged  beam  measurements  on  NOj  are  subject  to  some  important  uncer¬ 
tainties  that  are  not  present  in  measurements  on  other  ions.  First,  the  NO3  beam 
may  have  been  contaminated  with  some  NOJ  ions.  In  order  to  merge  NO*  and  NOJ 
at  low  relative  energy,  it  was  necessary  to  bring  the  negative  ion  beam  into  the 
merging  magnet  near  its  centerline.  This  reduced  the  total  mass  resolution  of  the 
system  to  a  point  where  NOj ,  which  is  formed  by  the  source  i.i  much  greater  abun¬ 
dance  than  NO3,  may  have  been  present  along  the  interaction  length.  Second, 
since  the  NO3  current  was  small,  the  resulting  neutral  signal  was  small,  and  it  is  pus- 
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sibic  thal  the  preamplifier  stage  u)  the  lock-in  amplifier  was  overloaded  by  the  non¬ 
coherent,  background  neutrals.  Third,  the  ion  NOj  exists  in  two  distinct  forms 
which  have  quite  different  electron  alTinities  (75) .  The  NOj  ions  used  in  the  merg¬ 
ed  beam  were  formed  in  a  mixture  of  NO  and  O2,  and  might  therefore  consist 
primarily  of  the  lower  electron  affinity  “peroxy”  form  (NO  •  Oj)  which  is  believed 
to  have  a  linear  configuration.  The  NOj  ions  of  Eisner  and  Hirsh  were  formed  in  an 
airlike  N2:  O2  mixture,  and  might  consist  primarily  of  the  “normal”  nitrate  form 
ofN03.  Tliese  two  forms  of  NOj  would  be  expected  to  have  different  neutraliza¬ 
tion  rates,  but  not  so  different  as  the  results  of  these  two  measurements.  A  fourth 
uncertainty,  which  stems  from  possible  vibrational  or  electronic  excitation  of  the 
ions,  is  common  to  all  measurements  using  molecular  beams,  and  will  be  discussed 
later. 

The  O2  +  NO3  reaction  has  also  been  investigated  using  the  merged  beam  appara¬ 
tus.  This  measurement  was  subject  to  the  same  uncertainties  as  the  NO'^  +  NOj 
measurement,  but  it  yielded  substantially  different  results.  The  mono-energetic  rate 
was  very  nearly  constant  from  0. 1 S  to  700  eV  at  a  value  of  1 .2,  and  the  extrapola¬ 
tion  yielded  a  thermal  energy  rate  of  ( 1.3  ±  0.4).  An  absorbing  sphere  calculation 
yielded  ( 1 .0  ±  0.2). 

7.4.  Disc  ussion  oj' molecular  ion  reactions 

The  thermal  reaction  rate  coefficients  obtained  for  all  systems  involving  mole¬ 
cular  ions  arc  summarized  in  table  2.  The  theoretical  results  are  all  from  absorbing 
sphere  calculations.  The  agreement  is  generally  within  a  factor  of  two  for  the  0~ 
and  O2  systems,  however  there  arc  larger  discrepancies  between  theory  and  experi¬ 
ment  for  NO^  +  0~  and  NO*  +  Oj  reactions  and  for  most  reactions  involving  NO2 
and  NO3.  In  general,  the  result  obtained  by  extrapolating  the  merged  beam  measure¬ 
ments  arc  larger  than  the  theoretical  predictions,  and  the  measurements  of  Eisner 
and  Hirsh  are  lower. 

It  is  possible  that  some  of  the  merged  beam  cross  sections  are  relatively  large  be¬ 
cause  of  internal  excitation  in  the  incident  ions.  The  ions  in  the  two  beams  are 
formed  in  duoplasmatron  ion  sources  and  probably  are  vibrationally  excited.  In  ad¬ 
dition,  if  metastable  electronic  states  exist  they  may  also  be  populated.  If  present, 
both  of  these  forms  of  excitation  would  tend  to  increase  the  neutralization  cross 
sections,  although  the  actual  magnitude  of  the  effect  is  unknown.  Electronic  ex¬ 
citation  generally  opens  up  more  possible  final  s.atcs,  and  the  general  effect  will 
be  to  increase  the  reaction  cross  sections.  Vibrational  excitation,  on  the  other  hand, 
has  a  more  limited  effect  since  it  will  increase  the  cross  sections  only  when  the  ef¬ 
fective  vibrational  overlap  with  the  final  states  is  increased  by  its  presence. 

From  the  observation  that  the  rate  for  NO*  and  any  negative  ion  were  found 
to  be  the  largest  for  that  negative  ion.  it  was  speculated  that  the  merged  beam  rates 
may  have  been  increased  by  the  presence  of  metastable  NO*  in  the  beam.  The  A^fl 
state  of  NO*  is  6  eV  above  the  ground  state  and  could  yield  much  larger  cross  sec¬ 
tions.  especially  in  the  case  of  NOj.  which  has  a  very  high  electron  detachment  en- 
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lablc  2 

Reaction  rate  coefficients  at  300°  K  -  inoicculai  ions.  o(300°  K)  -  1 0”’  cm*/sec 


System 

Experimental^ 

Theoretical  |42| 

llj  +  D~ 

4.7  ±  1.5 

8.5  1  2.1 

Nt  +0~ 

- 

2.0  1  0.6 

NO*  +  0“ 

4.9  t  2.0 

1 .9  1  0.6 

Oj  +0 

1.0  10.4 

1.9  1  0.5 

Nj 

1.6  10.5 

2.5  1  0.8 

NO*  +  Oj 

5.8  1  1.0 

2.4  1  0.8 

0*2  +07 

4.2  *  1.3 

2.4  1  0.8 

1.0  t0.1|74| 

Nt  +  NO7 

1.3  1O.5 

1.3  1  0.3 

NO*  +  NO7 

5  1  1  1.5 

1.2  1  0.3 

2.1  i0.6|27| 

1  75  1  0.6  |32| 

ot  +  N07 

4.1  1  1.3 

1.2  1  0.3 

NO*  +  NOJ 

8.1  i2.3 

1  1  1  0.3 

0.34  1  0  12  1 32) 

0*2  +  NO7 

1.3  1O.4 

1.0  1  0.2 

“  Experimental  results  ate  from  SRI  metpcd  beam  measurcmcnis  unless  otherwise  noted. 


ergy.  However,  nicasuremcitts  at  SRI  [73|  have  shown  lliat  less  than  lO'/f  of  the  NO* 
was  electronically  excited  in  the  merged  beam  experiments.  It  is  inconceivable  that 
such  a  small  excited  fraction  could  increase  the  rates  by  enougli  to  account  for  the 
observed  differences. 

For  NO3,  it  seents  likely,  as  is  mentioned  in  the  previous  section,  that  the  merged 
beam  measurements  were  made  using  the  low  electron  affinity  NO  •  0-,  form  of  the 
ion.  However,  it  is  doubtful  that  any  such  e.xcitation  can  account  for  more  than  a 
small  fraction  of  the  observed  differences.  For  example,  in  the  absorbing  sphere 
model,  decreasing  the  electron  attachment  energy  from  3.6  to  2.5  eV  increases 
the  calculated  thermal  energy  rate  by  only  \QI/i.  For  NOi.  a  decrease  from  2.3  to 
1.8  cV  increases  the  rate  by  only  1 57t. 

The  extrapolation  technique  dianissed  earlier  could  lead  to  larger  errors  than  ex¬ 
pected.  However,  the  Oi  +  NO^  and  NO*  +  NOT  rates  have  already  reached 
2.5  X  10"^  cm^/sec  at  0.1 5  eV  and  are  clearly  increasing  rapidly  at  the  low  energy 
limit  of  the  merged  beam  measurements.  It  ishiglily  unlikely  that  the  extrapolation 
technique  could  lead  to  errors  larger  than  a  factor  of  two  in  the  thennal  energy  reac¬ 
tion  rates. 

Tims  there  are  unresolved  discrepancies  between  experiments  and  theory  for 
rates  involving  NO2  and  NO3.  By  far  the  most  serious  discrepancy  is  l\)r  NO*  +  NO3. 
Because  of  the  importance  of  this  reaction  additional  investigations  are  planned  at 

SRI. 
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7,5,  Hydrated  ions 

Although  it  is  now  known  that  under  normal  quiescent  conditions  the  ions  in 
the  D-region  are  predominantly  in  hydrated  forms  as  discussed  brielly  in  section 
I ,  no  laboratory  measurements  have  yet  been  made  on  any  of  these  species.  Early 
efforts  at  SRI  were  thwarted  by  the  inability  to  produce  the  hydrated  species  in 
the  duoplasmalron  ion  sources.  The  ions  are  formed  in  three-body  reactions  such 
as  (7),  and  the  conditions  in  conventional  duoplasmatrons  apparently  either  are 
too  hot  or  too  violent  to  allow  production  in  sufficient  quantities  to  be  useful.  Re¬ 
cent  attempts  at  SRI,  using  a  duodehcairon  source  [76]  have  been  more  successful, 
and  there  now  is  a  good  chance  to  perform  measurements  on  at  least  some  of  the 
simpler  hydrates. 

As  water  molecules  are  attached  to  the  negative  ions,  several  effects  tend  to  slow 
down  the  neutralization  rate.  Firstly,  both  of  the  parent  ions  will  be  somewhat 
sliielded  from  each  other  by  the  water  molecules,  so  tliat  the  electron  transfer  pro¬ 
bability  is  reduced  at  each  moment  during  the  collision;  secondly,  the  energy  re¬ 
quired  to  remove  the  electron  from  the  negative  ion  is  increased  by  about  1.7  cV 
per  H2O  molecule,  so  that  the  number  of  neutral  state  crossings  available  to  the  re¬ 
action  is  decreased,  further  reducing  the  reaction  cross  section;  thirdly,  the  increas¬ 
ed  mass  reduces  the  average  velocity  at  each  temperature  and  tlie  rate  a  =  is  re¬ 
duced  even  if  the  reaction  cross  section  is  unchanged. 

On  the  other  hand,  another  factor  will  tend  to  increase  the  rate.  As  these  large 
molecules,  which  are  usually  polar  or  polarizable,  come  together  under  the  inlluence 
of  the  Coulomb  force,  internal  vibrational  and  rotational  modes  can  be  excited  by 
the  action  of  the  electric  field  on  the  moving  and  rotating  dipole.  If  the  induced 
e.xcitation  energy  in  these  modes  exceeds  AT.  the  ions  will  be  trapped  in  a  stable 
orbit,  and  many  more  curve  crossing  can  occur  to  allow  neutralization  before  the 
internal  excitation  can  be  removed  and  the  participles  fly  apart.  Or,  additional  in¬ 
ternal  modes  may  be  excited  to  bind  the  molecules  in  even  tighter  orbits.  For  very 
large  hydrated  ions  where  the  number  of  attached  water  molecules  exceeds  about 
5  or  more,  the  chance  for  orbital  capture  is  very  high  but  the  probability  of  elec¬ 
tron  transfer  may  be  very  low,  so  that  the  ions  many  actually  recombine  into  a 
large  neutral  water  cluster  bound  together  by  the  Coulomb  attaction  between  the 
core  ions.  Some  water  molecules  would  "boil”  off  during  the  process,  but  the  ions 
would  still  be  isolated  by  the  remaining  water  molecules. 

These  effects  arc  now  under  theoretical  study  at  SRI.  Huestis,  Smith  and  Ben¬ 
son  have  shown,  in  preliminary  results,  that  induced  rotational  excitation  can  af¬ 
fect  the  rates  even  for  a  simple  pair  such  as  O  +  NO'^  [  77] .  These  effects  have 
not  been  included  in  any  of  the  Landau-Zener  or  absorbing  sphere  calculations, 
and  may  account  for  some  of  the  tendency  for  tlie  theory  to  underestimate  the 
reaction  rates.  The  influence  of  internal  mode  excitation  should  be  considerably 
greater  for  the  hydrated  ions  than  for  the  simpler  ions  that  have  been  measured  so 
far,  thus  the  results  of  the  present  experiments  and  theoretical  efforts  should  be 
very  interesting. 
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8.  Concluding  remarks 

A  great  deal  of  progress  in  understanding  the  ion  ion  mutual  neutralization 
reaction  has  been  made  in  the  last  five  years.  The  merged  beam  experiments  have 
yielded  data  on  a  large  number  of  reactions  between  atomic  and  simple  molecular 
ions,  over  a  large  energy  range,  providing  a  very  useful  basis  for  the  testing  of  theo¬ 
retical  calculations,  as  well  as  data  from  which  could  be  obtained  reaction  rates 
that  were  previously  unknown.  These  rates  were  primarily  sought  for  use  in  iono¬ 
spheric  modeling,  but  they  are  also  important  for  the  understanding  of  flame  and 
gas  laser  phenomena.  A  theory  based  on  the  Landau  -Zener  approach,  using  semi- 
empirical  interaction  matrix  elements,  has  been  found  to  be  reasonably  good  for 
the  atomic  species,  although  it  cannot  predict  detailed  structure  such  as  found  ex¬ 
perimentally  in  the  H  case.  In  order  to  handle  the  more  complicated  prob¬ 
lem  of  molecular  ion  reactions,  the  absorbing  sphere  model  was  developed  and 
found  to  be  useful  for  estimating  unmeasured  reaction  rates  betwe  .-n  simple  mole¬ 
cules.  Thus  the  reaction  is  now  reasonably  well  understood  for  atomic  and  simple 
molecular  ions,  as  far  as  the  total  cross  section  or  reaction  rates  are  concerned. 

The  major  remaining  problem  from  the  practical  standpoint  of  the  reaction 
rates  needed  for  ionospheric  applications  is  to  determine  the  cross  sections  and 
rates  for  hydrated  ion  reactions.  Both  theoretical  and  experimental  efforts  are 
under  way  at  SRI,  aimed  toward  this  problem. The  theory  has  already  pointed  out 
the  importance  of  the  excitation  of  internal  energy  modes  during  the  collision, 
which  was  not  taken  into  account  in  the  treatment  of  the  simpler  molecules. 

Another  outstanding  problem  is  the  determination  of  the  final  states  and  pro¬ 
ducts  i)f  the  reactions.  Except  for  the  optical  measurements  on  Na"*^  +  0  ,  there 
have  been  no  attempts  to  study  this  aspect  of  ion  ion  neutralization.  Severe  ex¬ 
perimental  difficulties  are  associated  with  the  identification  of  the  neutral  product 
species  for  molecular  systems,  where  dissociation  is  often  likely  to  occur;  conse¬ 
quently  this  remains  as  a  major  unknown  in  atmospheric  modeling.  Because  of  the 
experimental  difficulties,  it  is  likely  that  the  best  information  on  the  reaction  final 
states  and  product  species  for  molecular  reactant  will  be  obtained  theoretically.  Pre¬ 
sent  theories  are  inadequate  for  such  detailed  predictions.  In  order  to  develop  ade- 
q.iate  theories,  more  detailed  experimental  data  are  required,  such  as  optical  mea¬ 
surements  of  final  states  of  reactions  between  atomic  ions  and  some  simple  mole¬ 
cular  systems. 
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